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AN  INVESTIGATION  OP  THE  ENERGY  CRITERION  IN  SNAP  BUCKLING 
PROBLEMS  WITH  APPLICATION  TO  CLAMPED  ARCHESf 

fey 

A,  Gjelsvik*  and  S.  R,  Bodner** 

ABSTRACT 

A  critical  analysis  Is  presented  of  the  significance  of 
the  energy  criterion  for  systems  that  exhibit  snap  buckling. 

The  results  of  energetic  analyses  are  compared  to  solutions  of 
the  differential  equations  for  the  nonlinear  behavior  of  a 
clamped  low  arch  subjected  to  a  central  concentrated  load.  The 
differential  equation  analyses  Include  consideration  of  transi¬ 
tional  non- symmetrical  buckling  modes.  The  snap  buckling  of 
low  arches  In  the  presence  of  high  temperature  creep  Is  also 
considered.  Experiments  are  described  on  the  deformation  of 
such  arches  under  static  and  creep  conditions. 


*  The  results  in  this  paper  were  obtained  in  the  course  of 
research  sponsored  by  the  Office  of  Naval  Research  under 
Contract  Nonr-£62(20)  with  Brown  University, 
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Principal  Symbols 


central  height  of  mechanical  model;  also  non- d Imena 1 onal 
amplitude  of  symmetrical  deflection  of  arch 

non-dimensional  amplitude  of  antl-symmetrlcal  deflection 
of  arch 

width  of  arch  cross  section 
spring  constant 
thickness  of  arch 

tangential  and  radial  displacements  at  middle  surface 
of  arch 

generalized  displacement 

cross  sectional  area;  also  deflection  amplitude 
deflection  amplitude 
Young* a  modulus 

force  In  spring  component  of  the  model 

total  energy  of  system 

HX/P5 

applied  load 

energy  load  (same  as  "Intermediate”  or  "dead  weight" 
buckling  load) 

lower  buckling  load 

upper  buckling  load 

FK/Et2fp 

length  of  arms  of  model 

radius  of  arch  centerline 

axial  force  in  the  arms  of  the  model 

angular  coordinate  for  arch  (measured  from  centerline) 
half  of  angle  subtended  by  arch 


Y  Imperfection  amplitude 

6  deflection  of  midpoint  of  model 

e  strain 

k  *  «/p 

x  change  of  curvature  of  arch 

X  *  p2R/t 

d  stress 
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Introduction. 

The  fundamental  cause  of  the  large  discrepancy  between  buck¬ 
ling  values  predicted  by  the  classical  eigenvalue  theory  and  exper¬ 
imental  results  for  spherical ,  and  cylindrical  shells  was  shown  by 
von  Karman  and  Tsien,  [1]  [2],  to  be  due  to  the  existence  of  post 
buckling  stable  equilibrium  states  characterized  by  loads  appreci¬ 
ably  lower  than  the  classical  buckling  load,  e.g.  Pig.  1,  They 
advanced  the  hypothesis  that  somehow  the  shell  could  jump  to  these 
equilibrium  states  at  loads  less  than  the  linear  buckling  load  so 
that  the  lowest  load  for  such  states,  point  2  of  Pig.  1,  should  be 
considered  the  minimum  load  at  which  buckling  was  possible.  The 
question  was  left  open  as  to  how  the  transition  from  the  initial 
state  to  the  buckled  state  was  to  take  place  since,  unlike  classical 
buckling  theory,  these  are  not  geometrically  adjacent  states.  They 
postulated,  however,  that  the  transition  may  be  connected  with 
geometrical  imperfections,  small  dynamic  disturbances  or  possibly 
with  unsymmetrical  deformations  not  included  in  the  analyses  [3]# 

The  values  von  Karman  and  Tsien  obtained  for  the  lower  buckling 
loads  in  [1]  and  [2]  were  in  surprisingly  good  agreement  with  avail¬ 
able  experimental  results, 

A  weakness  in  the  lower  buckling  load  as  a  realistic  buckling 
criterion  was  pointed  out  by  Friedrichs  [I4.] .  He  showed  that  at  the 
lower  buckling  load  the  system  of  the  shell  and  the  loading  is  at 
a  higher  energy  level  in  the  buckled  than  in  the  unbuckled  state. 
Since  no  energy  is  transferred  to  the  system  during  the  jump,  but 
3ome  may  be  lost,  it  can  only  be  expected  to  Jump  from  a  state  of 
a  higher  to  a  state  of  lower  energy  level.  If  no  energy  is  lost 
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the  minimum  load  at  which  a  jump  can  take  place  is  then  the  load  at 
which  the  total  energy  in  the  buckled  and  unbuckled  states  are  equal. 
Friedrichs  called  this  the  Intermediate  buckling  load.  He  did  not 
try  to  explain  why  or  how  the  jump  could  take  place.  Friedrichs 
also  analyzed  the  stability  of  the  complete  sphere  using  a  boundary 
layer  treatment  and  relaxed  some  unnecessarily  restrictive  assump¬ 
tions  employed  by  von  Karman  and  Tsien.  Friedrichs  found  that  both 
the  minimum  and  the  intermediate  buckling  loads  were  zero  in  the 
boundary  layer  development  for  a  complete  sphere. 

Tsien  [5,6]  attempted  to  explain  the  mechanism  of  snap  buck¬ 
ling  by  introducing  additional  energy  considerations.  His  energy 
criterion  stated  that  under  average  laboratory  or  actual  service 
conditions  the  most  probably  equilibrium  state  is  that  with  the 
lowest  possible  energy  level.  In  calculating  the  energy  levels  be¬ 
fore  and  after  buckling  the  elasticity  of  the  loading  system  must 
be  taken  into  account.  Different  buckling  loads  are  therefore  pre¬ 
dicted  for  the  sfv'ne  shell  with  different  loading  systems  by  this 
theory.  For  dead  loading  the  load  necessary  to  cause  instability 
is  the  same  as  Friedrich's  intermediate  buckling  load.  Tsien* s 
principle,  however,  adds  the  assumption  that  there  always  are  dis¬ 
turbances  of  sufficient  magnitude  present  to  overcome  the  stability 
of  the  unbuckled  state.  However,  if  disturbances  of  such  a  magni¬ 
tude  are  admitted  there  appears  to  be  no  reason  why  the  system 
could  not  jump  from  a  state  of  lower  to  a  state  of  higher  energy, 
provided  a  stable  buckled  state  exists  at  that  load.  The  lower 
buckling  load  is  the  lowest  load  for  which  a  stable  buckled  state 
exists,  and  therefore  the  lowest  load  for  which  the  system  could 
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a 

remain  in/buckled  position  If  sufficiently  disturbed.  Below  this 
load  only  oscillations  about  the  unbuckled  state  can  take  place. 

The  energy  load,as  defined  In  this  manner,  therefore  loses  Its 
original  meaning, 

Tslen  Investigated  the  stability  of  complete  spheres  under 
hydrostatic  pressure  and  cylindrical  shells  under  axial  compression, 
and  obtained  fairly  good  agreement  with  experimental  results.  Al¬ 
though  Tslen* s  results  and  more  recent  work  using  the  energy  cri¬ 
terion  by  Lo,  Crate  and  Schwarz  [7]  gave  good  agreement  with  exper¬ 
iments,  there  are  strong  objections  to  Its  validity  and  therefore 
to  Its  acceptability.  These  objections  are  well  summarized  in  the 
survey  paper  by  Fung  and  Sechler  [8]  and  are  based  on  the  fuzziness 
of  the  theoretical  foundation  of  the  energy  criterion  and  especially 
the  fact  that  it  requires  a  knowledge  of  the  geometrically  removed 
post  buckling  state  to  determine  the  onset  of  buckling.  In  addi¬ 
tion,  the  t e st s  of  Kaplan  and  Fung  on  shells  [9]  showed  little 
difference  for  different  loading  conditions  In  opposition  to  the 
theory.  Tests  on  simply  supported  shallow  arches  [10]  also  agreed 
well  with  the  upper  buckling  load  values  which  were  appreciably 
different  from  the  loads  predicted  by  the  energy  criterion. 

At  the  present  time  the  trend  Is  away  from  Intuitive  ap¬ 
proaches  such  as  the  energy  criterion  and  attempts  are  made  to  solve  the 
governing  differential  equations  by  the  use  of  high  speed  digital 
computers,  e.g,  Keller  and  Reiss  [11],  Budlansky  [12]  and  Thurston 
[13],  There  Is  little  doubt  that  this  approach  will  eventually 
give  most  of  the  required  answers  since  for  completely  static  load¬ 
ing  the  buckling  load  must  be  a  peak  point  of  the  appropriate  load 


deflection  curve.  Such  curves  are,  however,  extremely  difficult 
to  obtain  since  they  should  incorporate  initial  deflections  of  the 
shell,  the  appropriate  stresses  and  deformations  that  exist  in  the 
pre-buckling  state  and,  above  all,  must  include  the  correct  mode 
into  which  the  shell  would  initially  buckle.  Aside  from  the  basic 
problems  associated  with  handling  nonlinear  problems,  additional 
complications  arise  when  the  important  initial  deflections  are  not 
necessarily  those  having  the  form  of  the  buckled  shape, and  when 
transitional  buckling  modes  appear  that  may  also  be  appreciably 
different  from  both  the  initial  and  the  final  shapes.  The  formula¬ 
tion  and  numerical  solution  of  nonlinear  problems  of  this  generality 
still  presents  formidable  obstacles. 

Because  of  the  difficulties  mentioned  above  in  contrast  to 
the  simplicities  offered  by  an  energy  approach,  the  authors  have 
felt  it  worthwhile  to  reexamine  the  theoretical  basis  of  the  energy 
criterion.  Its  relative  simplicity  makes  it  particularly  attractive 
for  problems  in  which  additional  complicating  features  such  as 
anisotropy,  inhomogeneity,  and  combined  loadings  appear.  In  the 
following  discussion  the  term  "energy  load"  or  "energy  criterion" 
will  be  restricted  to  mean  I,Yiedrich|s  intermediate  load  or,  alter¬ 
natively,  Tslen*s  "dead  weight"  energy  load. 

The  first  part  of  this  paper  is  devoted  to  an  analysis  of 
the  significance  of  the  energy  load  for  systems  having  force  de¬ 
flection  characteristics  of  the  form  shown  in  Fig,  1,  The  general 
disousslbn  is  preceded  by  a  complete  analysis  of  a  simple  mechani¬ 
cal  model  which  has  such  a  load  defleotion  curve.  The  discussion 
Itself  is  completely  general  and  is  not  restricted  to  arch  and  shell 


structures  -  the  arguments  may  be  relevant  to  other  systems  which 
experience  similar  behavior. 

The  second  part  of  this  paper  Is  devoted  to  comparing  the 
results  of  the  energy  approach  with  those  of  more  thorough  theore¬ 
tical  analyses  of  low  arches.  A  subsequent  paper  will  Investigate 
the  energy  criterion  in  relation  to  the  behavior  of  clamped 

spherical  caps  and  complete  spherical  shells. 

A  fairly  complete  theoretical  and  experimental  analysis  of 
clamped  shallow  arches  Is  performed.  The  theoretical  work  includes 
consideration  of  a  number  of  symmetrical  deflection  modes  as  well 
as  transitional  non- symmetrical  buckling.  The  experiments  on  arohes 
produced  buckling  load  values  and  complete  load  deflection  curves 
which  Included  the  unstable  transition  region.  The  agreement  be¬ 
tween  theory  and  experiment  was  fairly  good. 

The  energy  load  computed  on  the  basis  of  Ideal  geometry  and 
axially  symmetrical  buckling  Is  shown  to  be  a  lower  bound  for  snap 
buckling  for  systems  influenced  by  non-symmetrical  transitional 
buckling  modes.  For  a  particular  class  of  systems,  the  energy  load 
Is  also  found  to  be  unaffected  by  certain  forms  of  initial  Imper¬ 
fections  and  to  be  a  lower  bound  for  snap  buckling  In  the  presence 
of  those  imperfections.  The  systems  and  imperfections  are  restrict¬ 
ed  to  those  for  which  the  load  deflection  curve  with  and  without 
the  inclusion  of  imperfections  is  anti- symmetrical  about  a  common 
straight  line  through  the  origin.  Although  only  few  structural 
systems  fall  into  this  class,  e.g.,  clamped  arches,  the  actual 
behavior  of  spherical  caps  and  presumably  other  systems  Is  close  to 
being  the  same.  The  energy  load  Is  found  to  be  relatively 
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Insensitive  to  arbitrarily  chosen  buckling  modes.  In  addition,  both 
the  theoretical  and  experimental  results  for  arches  ar.d  spherical 
caps  show  that  the  difference  between  the  energy  load  and  the  upper 
buckling  load  diminishes  the  more  the  system  is  influenced  by 
non- symmetrical  transitional  modes  and  by  destabilizing  initial 
imperfections. 

Model  for  Snap  Buckling, 

To  gain  a  better  understanding  of  some  of  the  concepts 
mentioned  in  the  introduction,  as  well  as  to  introduce  some  new 
ones,  a  simple  mechanical  model  will  first  be  investigated.  The 
bar  spring  system  shown  in  Fig,  2  has  many  important  features  in 
common  with  a  low  arch  and  a  shallow  cap.  The  model  is  considered 
as  hinged  in  the  middle  and  at  the  ends,  and  subjected  to  a  single 
load  P  at  the  apex.  The  discussion  will  be  limited  to  a  low  frame, 
that  is  one  for  which 


(a/L)«  1  (1) 

where  a  is  the  unstrained  height  and  2L  is  the  width  of  the  frame. 
Prom  the  geometry  the  equilibrium  equation  is  immediately  obtained: 

P  =  kb  +  2T(a-6)/L  (2) 

where  6  is  the  vertical  displacement  of  the  apex,  k  the  spring 
stiffness  and  T  the  compressive  force  in  the  arms.  Also  from 
geometry  the  strain  in  the  arms  is 

e  =  6(6-2a)/2L2 


(3) 
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Prom  Hook’s  law  the  force  in  the  arms  is  therefore 

T  *  -  (AE/2L2)6(6-  2a)  (i*) 

where  AE  is  the  stiffness  of  the  arms.  Substituting  for  T  in  the 
equilibrium  equation  (2)  yields  the  load  deflection  relation  of 
the  model 

P  =  k6  +  (AE/L3)6(6-a)  (&-2a)  (5) 


This  is  plotted  in  Fig.  3. 

The  curve  can  be  considered  as  composed  of  two  parts,  a 
straight  line  caused  by  the  spring  stiffness  k,  and  a  curve  anti- 
symmetrical  about  6  =  a.  According  to  the  classical  criterion 
instability  occurs  when  the  first  maximum  is  reached  or  at  the 
point  A  in  Pig,  3,  This  is  usually  called  a  limit  or  critical 
point.  In  contrast  to  most  buckling  problems  instability  is  not 
caused  by  a  bifurcation  of  the  equilibrium  state  but  by  the  load 
reaching  a  maximum.  The  lower  buckling  load  of  von  Karman  and 
Tsien  would  correspond  to  the  minimum  point  or  point  B  in  Pig,  3* 
Both  these  points  are  obtained  by  setting 


dP  .. 


0 


(6) 


By  substitution  from  (5),  this  yields 

6  *  a  +  [a2/3)-(kL3/3AE)]* 


(7) 


The  plus  sign  la  clearly  associated  with  the  lower  buckling  load 
and  the  negative  with  the  limit  point  or  the  upper  buckling  load. 
These  two  loads  will  coincide  in  a  point  of  inflection  with  a 
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horizontal  tangent  at  6  =  a  when  the  quantity  under  the  root  ia 
zero.  If  the  root  is  imaginary  the  load  deflection  curve  will  be 
continually  rising  and  always  stable,  For  instability  to  be  possi¬ 
ble  therefore 

(a2AE/kL3)  >1  (8) 

The  intermediate  buckling  load  or  the  energy  load  as  defined 
by  Friodriohs,  or  by  Tslen  for  a  dead  load,  can  be  determined  most 
easily  from  the  geometry  of  the  load  deflection  curve.  At  the 
energy  load  the  total  energy  in  the  buckled  and  unbuckled  states 
are  equal,  or  the  change  in  energy  from  one  to  the  other  is  zero. 

In  Fig,  4-  the  decrease  in  potential  energy  of  the  load  in  going 
from  the  unbuckled  to  the  buckled  state  is  given  by  the  rectangle 
ABCD,  The  Increase  in  strain  energy  of  the  spring  and  the  arms  is 
given  by  the  area  under  the  load  deflection  curve  between  B  and  C, 
The  change  in  total  energy  will  then  be  given  by  the  difference 
between  these  two  areas.  At  the  energy  load,  the  two  areas  must 
be  equal.  From  the  symmetry  of  the  curve  it  is  seen  that  this 
occurs  when  the  line  BC  passes  through  the  crossing  point  of  the 
straight  line  due  to  the  spring  and  the  curve  due  to  the  effect  of 
the  arms.  The  energy  load  is  therefore 

Pe  =  ka  (9) 

This  ia  a  function  of  the  spring  stiffness  and  the  geometry  but  is 
independent  of  the  properties  of  the  arms.  This  is  quite  remarka¬ 
ble  and  shows  that  imperfections  thathave  the  effect  of  reducing 
the  stiffness  of  the  arms  only  will  leave  the  energy  load  unchanged. 
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The  energy  load  will  cease  to  exist  when  the  curve  is  reduced  to 
one  with  a  point  of  inflection  with  a  horizontal  tangent.  The 
upper,  lower  and  energy  load  would  then  be  equal.  If  imperfections 
of  the  arms  are  sufficiently  large,  the  energy  load  may  cease  to 
exist,  but  in  that  case  instability  would  not  occur  either.  Con¬ 
sideration  of  imperfections  of  the  arms  only  would  result  in  the 
enerry  load  bei.ne.  a  lo^nr  bound  on  the  loads  that  can  cause 
instability.  This  is  true  regardless  of  the  magnitude  of  the  im¬ 
perfections.  It  may  be  objected  that  the  height  "a"  of  the  model 
enters  into  the  expression  for  the  energy  load.  Variations  in  "a” 
however,  should  not  be  considered  to  be  imperfections  in  the  usual 
sense  of  buckling  problems  but  as  changes  in  overall  geometry. 

The  same  could  be  said  about  the  length  of  a  strut  in  ordinary 
Euler  buckling^ i.e.  imperfections  are  considered  to  be  deviations 
from  straightness  and  not  changes  in  the  length.  In  the  general 
class  of  nonlinear  problems  under  consideration,  however,  more 
attention  must  be  given  to  the  definition  of  geometric  "imperfec¬ 
tions". 


As  an  example  of  an  imperfection,  the  arms  can  be  considered 
to  have  a  sinusoidal  centerline  with  amplitude  a  .  The  shortening 
e  of  one  of  the  arms  subjected  to  the  compressive  force  T  is  then: 


e 


2  ip  i 

^"^-(Ai/l2)1 


2 

rc 

T* 


TL 

aE 


(10) 


This  must  be  equal  to  the  shortening  of  the  arms  as  caused  by  the 
lowering  of  the  apex  of  the  model. 

Then  from  Eq.  (3), 


(11) 
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Prom  the  equilibrium  equation  (2)  the  value  of  T  can  be  substituted 
to  give  the  load  deflection  curve 


v  252T  (P-kM/U-6)  I2 

^t(P-k4)/(a-6)1.2^EX/L3j 


( P-k6 ) L2 
2AE(  a-6  )** 


6(2a-6) 

“TE 


=  0 


(12) 


This  is  plotted  in  Pig,  5  for  different  values  of  y  ,  It  can  be 
seen  that  the  upper  buckling  load  is  reduced  by  y  increasing  where¬ 
as  the  energy  load  remains  unchanged. 

If  the  possibility  of  the  arms  buckling  is  admitted  a  dif¬ 
ferent  mode  of  instability  is  added.  Let 


Te  =  n2El/L2 


(13) 


be  the  Euler  load  of  the  arms,  Por 

T  <Te 

the  load  deflection  curve  is  still  given  by  (5),  but  when  T  =  Tg 
it  i 3  given  by  substituting  for  T  from  (13)  in  the  equilibrium 
equation  (2),  This  yields 

P  =  k&  +  (a-6)  ( II4.) 

which  is  a  straight  line  passing  through  P  =  ka  and  6  =  a  as  shown 
in  Pig,  6,  Instability  will  now  occur  when  the  force  in  the  arms 
becomes  Tg  or  at  the  point  B  in  the  figure  provided  the  slope  of 
the  line  BC  is  negative  or  zero.  If  the  slope  is  positive  the 
load  deflection  curve  will  be  continually  rising  and  stable  for 
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all  loads,  B’rom  Fig,  6  it  is  seen  that  the  slope  becomes  positive 
when  the  column  buckling  load  is  less  than  the  energy 
load  P  ,  The  energy  load  is  therefore  again  the  dividing  point 
between  stability  and  instability  if  the  possibility  of  transi¬ 
tional  buckling  modes  is  admitted.  The  energy  load  is  therefore 
also  in  this  sense  a  lower  bound  on  the  snap  buckling  load  of  the 
complete  system. 

From  the  discussion  of  the  model  it  is  seen  that  the  energy 
load  calculated  on  the  basis  of  perfect  geometry  and  the  basio 
deformation  mode  is  a  lower  bound  on  the  snap-buckling  load  in 
the  presence  of  imperfections  or  additional  transitional  modes.  It 
is  of  interest  to  see  if  this  interpretation  of  the  energy  load  is 
true  only  for  this  particular  model  or  if  it  is  true  for  more 
general  classes  of  systems. 

Influence  of  Transitional  Buckling  and  Imperfections  in  Systems 
fhat  fixhlbH  Snap  Suckling,  ’ 

For  systems  having  load  deflection  curves  similar  to  Fig.l 
in  the  absence  of  both  transitional  buckling  modes  and  initial  im¬ 
perfections,  the  problem  is  to  determine  the  changes  those  factors 
could  have  on  the  basic  curve  and  tne  relevancy,  if  any,  of  the 
energy  criterion.  The  nature  of  admissible  load  deflection  curves 
for  transitional  buckling  modes  will  first  be  examined.  Transi¬ 
tional  buckling  refers  to  the  case  where  the  mode  of  the  inter¬ 
mediate  unstable  equilibrium  states  is  essentially  different  from 
both  the  initial  state  and  the  final  buckled  state.  Examples  of 
transitional  buckling  are  column  buckling  of  the  model  previously 
discussed  and  anti-symmetrical  buckling  of  simply  supported  arches 
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which  was  investigated  by  Marguerre  [141  and  is  discussed  in  [3]. 

The  basic  system  under  discussion  will  be  assumed  to  have 
only  one  pair  of  stable  initial  and  final  equilibrium  states  in 
the  load  range  <  P  _<  as  indicated  by  the  curves  OA  and  BC 
respectively  in  Fig.  7*  Load  deflection  curves  that  cross  over 
themselves  will  be  ruled  out.  The  overall  (or  average)  behavior 
of  most  structural  elements  that  exhibit  snap  buckling  would  have 
equilibrium  characteristics  similar  to  Fig.  7.  The  relevant 
question  for  transitional  buckling  is  the  nature  of  possible  un¬ 
stable  equilibrium  paths  connecting  the  two  stable  branches.  Alter¬ 
natively,  the  problem  is  the  determination  of  admissible  paths 
connecting  the  two  stable  states  that  do  not  require  any  input  or 
release  of  kinetic  energy. 

In  discussing  snap  buckling  of  a  system  with  many  degrees 
of  freedom  It  Is  necessary  to  introduce  some  characteristic  dis¬ 
placement  function  in  order  to  describe  the  geometric  dependence. 
The  function  that  appears  to  be  an  appropriate  and  consistent  in¬ 
dicator  of  the  gross  behavior  is  the  one  that  appears  in  the  ex¬ 
pression  for  the  work  done  by  the  external  force  system  (the 
generalized  displacement).  For  the  case  of  a  concentrated  force 
this  would  be  the  actual  displacement  at  the  point  of  application 
and  in  the  direction  of  the  applied  force.  For  hydrostatic  pres¬ 
sure  it  would  be  the  change  in  volume.  The  integrated  average 
displacement  used  as  the  geometric  reference  by  Budiansky  [12]  is, 
to  first  order  terms,  equivalent  to  a  volume  change.  The  general¬ 
ized  displacement  from  the  work  expression  will  be  used  as  the 
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geometric  reference  parameter  In  the  following  general  discussions 
and  In  the  particular  examples. 

The  elastic  strain  energy  U  of  the  system  Is  assumed  to  be 
completely  defined  In  the  unbuckled  and  buckled  stable  equilibrium 
states, The  work  done  by  the  external  force  system  P  during  a  small 
change  In  the  generalized  displacement  v  Is  Pdv,  In  the  absence 
of  klnetlo  energy  this  must  be  equal  to  the  change  of  c'ure.in  energy 
during  the  same  process, 

dU  ■  Pdv  (15) 

If  the  total  potential  energy  H  defined  as 

H  =  U  -  Pv 


Is  introduced,  then 


dH  =  dU  -  Pdv  -  vdP 


and  from  (15) 


dH  =  -  vdP 


(16) 


Prom  the  definition  of  the  total  energy  It  is  seen  that  It  Is 
also  completely  defined  In  the  unbuckled  and  buckled  stable  equi¬ 
librium  states.  Integrating  (16)  along  a  path  connecting  two 
such  states  a  and  b  gives 

b 

Hb  -  Ha  =  -  i  vdP 

a 


(17) 
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where  the  integral  depends  upon  the  path.  This  integral  can  be 
interpretated  graphically  as  the  area  between  the  path  and  the  P 
axis.  This  area  must  therefore  be  the  same  for  all  paths  between 
the  states  a  and  b  if  no  kinetic  energy  is  involved.  At  the 

energy  load,  Pft,  in  particular,  =  Hg,  and  setting  a  =  1,  b  =  2 
to  correspond  to  the  points  in  Pig,  7, 

2 

JvdP  =  0  (18) 

1 

and  the  area  is  therefore  zero.  This  is  the  same  criterion  that 
was  used  to  determine  the  energy  load  for  the  model.  Some  equili¬ 
brium  curves  that  would  satisfy  (18)  are  shown  in  Pig,  7  by  the 
dashed  lines.  Portions  of  these  curves,  except  the  straight  line, 
would  be  stable  equilibrium  states  aside  from  OA  and  BC  and  must 
therefore  be  excluded.  The  straight  horizontal  line  is  therefore 
the  only  possible  completely  unstable  or  neutral  equilibrium  path 
connecting  the  unbuckled  and  buckled  states  at  the  energy  load. 
Starting  from  above  the  energy  load,  possible  unstable 
equilibrium  must  be  similar  to  the  curve  3-4  l-1*  Pig*  7*  Since 
the  portion  of  the  curve  1-3-5-4-2  can  be  considered  as  a  path 
connecting  1  and  2  for  which  (18)  holds,  it  follows  that  the  area 
enclosed  by  1-3-5  must  be  equal  to (minus)  the  area  enclosed  by 
5-4-2,  The  path  3-4  must  also  be  monotonically  decreasing  to 
satisfy  the  condition  that  it  be  a  locus  of  unstable  or  neutral 
equilibrium  states.  Starting  from  points  below  the  energy  load, 
no  path  satisfying  the  requirements  can  be  found.  The  energy  load 
is  therefore  the  minimum  load  from  which  unstable  equilibrium 
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paths  could  start,  i,e,  It  Is  the  minimum  load  for  which  a  bifurca¬ 
tion  of  equilibrium  due  to  the  occurrence  of  other  modes  would 
lead  to  snap  buckling  of  the  system. 

Some  other  general  statements  follow  almost  directly.  For 
a  system  with  a  load  deflection  curve  of  the  type  shown  in  Fig,  1 
the  energy  load  must  lie  between  the  upper  and  lower  buckling  loads 
since  the  two  shaded  areas  must  be  equal.  The  three  critical  loads 
are  equal  when  the  load  deflection  curve  has  a  point  of  inflection 
with  a  horizontal  slope.  It  is  also  seen  that  for  the  same  load 
the  total  energy  H  along  BC  (Fig,  7)  is  less  than  H  along  OA  for 
P  >  P ,  Since  BC  is  the  locus  of  stable  equilibrium  states,  the 

w 

potential  energy  in  adjacent  non-equilibrium  states  (at  the  same 

load)  must  be  greater  than  for  points  on  BC,  The  energy  load  P 

© 

can  therefore  be  defined  as  the  minimum  load  for  which  the  total 
potential  energies  are  equal  in  the  unbuckled  state  (along  OA)  and 
in  any  arbitrary  buckled  state  whether  or  not  it  is  an  equilibrium 
state.  The  buckled  state  that  gives  the  true  minimum  must,  of 
course,  be  the  equilibrium  state  2,  This  variational  formulation 
for  calculating  the  energy  load  may  be  useful  since  it  does  not 
require  the  prior  determination  of  the  buckled  equilibrium  states 
BC,  The  energy  load  determined  by  such  a  variational  procedure 
would  be  an  upper  bound  if  the  initial  state  along  OA  is  exact  and 
the  buckled  state  is  approximated,  and  would  be  a  lower  bound  on 
the  real  energy  load  if  the  reverse  were  true. 

The  role  of  imperfections  in  nonlinear  systems  is  a  much 
more  difficult  problem.  Part  of  the  difficulty  is  due  to  the 
vagueness  of  any  definition  of  "imperfections"  other  than  any 
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arbitrary  deviation  from  the  pre-assumed  "Ideal”  geometry  that 
satisfies  the  geometric  boundary  conditions.  It  Is  not  surprising 
that  "Imperfections"  can  have  vastly  different  effects  according 
to  their  manner  of  Introduction  In  a  specific  problem. 

In  linear  buckling  problems,  imperfections  are  generally 
taken  to  mean  the  deviations  from  the  idealized  geometry  that  in¬ 
fluence  the  destabilizing  (or  stabilizing)  role  of  the  membrane 
stresses.  This  same  concept  Is  usually  used  in  nonlinear  problems 
In  the  sense  that  the  terms  due  to  the  Imperfections  appear  only 
in  expressions  for  the  membrane  energy  or  In  conjunction  with  mem¬ 
brane  stresses  In  the  differential  equations.  The  Imperfections 
do,  however,  influence  the  complete  problem  since  the  membrane 
and  bending  stresses  and  the  displacements  are  all  coupled  in  a 
nonlinear  problem. 

The  imperfections  considered  in  the  analysis  of  the  model 
were  deviations  from  straightness  of  the  arms  and  therefore  had 
the  same  mode  as  that  corresponding  to  transitional  buckling. 

The  energy  load  was  found  to  be  unchanged  and  to  be  a  lower  bound 
for  snap  buckling  in  the  presence  of  these  imperfections.  Another 
kind  of  imperfection  in  the  model  problem  would  be  a  change  in  the 
height  of  the  apex  which  would  correspond  to  a  change  In  the  basic 
geometry  and  therefore  in  the  energy  load.  It  is  convenient  to 
designate  these  two  basically  different  types  of  imperfections  as 
A  and  B  respectively  for  reference  in  discussing  more  complicated 
systems. 

The  general  conditions  req>ired  for  imperfections  to  be  of 
type  A  are  of  particular  Interest  since  the  energy  criterion  then 


17 


has  a  definite  physical  interpretation.  We  shall  consider  systems 
for  which  the  resistance  to  deformation  can  be  subdivided  into 
"bending"  and  "membrane"  contributions  and  where  the  membrane 
stresses  are  symmetrical  in  the  stress  -  generalized  displacement 
plane  with  respect  to  a  particular  displacement  value.  This  char¬ 
acteristic  is  shown  in  Fig,  8  where  the  stress  deflection  rela¬ 
tionships  are  assumed  to  be  symmetrical  about  the  vertical  line  A- A, 
The  strain  energy  of  the  membrane  stresses  is  a  quadratic  function 
of  the  stresses  and  would  therefore  have  the  general  form  shown  in 
Fig.  9.  It  would  also  be  symmetric  about  A-A,  The  membrane  stress 
contribution  to  the  load  deflection  curve  would  then  be  anti-sym¬ 
metric  about  A-A,  since  that  curve  is  obtained  by  differentiating 
the  potential  energy  with  respect  to  the  generalized  deflection  v  , 
In  nonlinear  structural  systems  the  terms  corresponding  to  the 
bending  resistance  are  generally  linear  in  the  displacements.  The 
complete  load  deflection  curve  would  then  have  the  same  form  as 
that  for  the  model.  Fig,  3« 

If  the  inclusion  of  geometric  imperfections  into  the  problem 
only  alters  the  magnitude  of  the  membrane  stresses  but  does  not 
change  the  bending  contribution  or  the  basic  symmetry,  then  the 
energy  load  must  remain  unchanged,  A  sufficient,  and  seemingly 
necessary,  condition  that  the  energy  load  remain  unchanged  in  the 
presence  of  geometric  imperfections  is,  therefore,  that  the  membrane 
stresses  for  the  perfect  system  are  all  symmetrical  about  a  common 
generalized  displacement  value,  and  that  the  only  effect  of  the  im¬ 
perfections  considered  is  to  change  the  magnitudes  of  those  membrane 
stresses  without  altering  the  symmetry,  A  further  consequence  of 
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these  oonditions  Is  that,  as  was  shown  In  the  case  of  the  model, 
the  energy  load  would  also  be  a  lower  bound  on  snap  buckling  in 
the  presenoe  of  the  imperfections.  The  imperfections  would  then 
be  of  type  A  discussed  previously. 

It  does  not  seem  possible,  however,  to  determine  a  priori 
whether  a  given  structural  system  and  arbitrarily  chosen  imperfec¬ 
tion  functions  will  behave  in  the  above  manner.  In  the  subsequent 
analyses  it  is  found  that  the  behavior  of  the  symmetrically  loaded 
clamped  arch  does  eonform  with  the  stated  conditions.  For  the 
clamped  aroh  anti- symmetrical  imperfections  are  found  to  be  of 
type  A  (the  transitional  buckling  mode  is  also  anti- symmetrical ) , 
while  symmetrical  imperfections  are  of  type  B.  The  clamped  spheri¬ 
cal  cap,  on  the  other  hand,  does  not  possess  the  required  symmetry 
properties  of  the  membrane  stresses,  so  that  it  cannot  admit  im¬ 
perfections  strictly  of  the  A  type.  The  symmetrical  imperfections 
considered  by  Budiansky  in  [12]  for  the  cap  appear  to  have  proper¬ 
ties  similar  to  those  of  type  B, 

Cylindrical  shells  also  do  not  appear  to  have  the  properties 
required  for  a  class  of  imperfections  to  be  strictly  of  type  A. 

It  is  interesting  to  note,  however,  that  the  results  of  [15,  16] 
approximate  the  expected  results  for  imperfections  of  type  A,  The 
upper  buckling  load  for  the  largest  imperfection  parameter  that 
shows  snap  buckling  is  reasonably  close  to  the  energy  load  of  the 
perfect  shell. 

Although  the  energy  criterion  does  not  maintain  any  definite 
physical  meaning  in  the  presence  of  imperfections  of  type  B,  which 
can  be  related  to  changes  in  the  overall  geometry,  it  can 
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serve  as  a  measure  of  the  gross  influence  of  these  imperfections 
on  the  nonlinear  behavior  of  the  system.  The  change  in  the  energy 
load  due  to  those  Imperfections  would  be  an  indication  of  the  sen¬ 
sitivity  of  the  system  to  the  geometry  changes.  This  viewpoint  is 
of  some  practical  interest  since  the  computation  of  the  energy 
load  is  considerably  simpler  than  the  complete  solution  of  the  non¬ 
linear  differential  equations.  This  would  be  particularly  true 
in  problems  where  the  mode  of  sensitive  imperfections  of  type  B 
are  essentially  different  from  either  the  transitional  or  final 
buckled  modes,  e,g,  axi-symmetrical  imperfections  in  cylinder 
buckling  problems. 

Analysis  of  a  Clamped  Low  Arch  Subjected  to  a  Central  Concentrated 
Load 

The  clamped  low  arch  was  chosen  as  a  specific  structural 
system  to  examine  the  general  results  of  the  preceding  section. 
Although  the  simply  supported  arch  has  been  thoroughly  investi¬ 
gated  by  Fung  and  Kaplan  [10]  and  by  Biezeno  and  Grammel  [17],  the 
clamped  condition  Initially  appeared  to  offer  more  basic  compli¬ 
cating  features.  In  nonlinear  problems  particularly,  a  change  in 
the  boundary  conditions  can  often  alter  the  complete  nature  of  the 
problem.  The  anticipated  difficulties  did  not,  in  fact,  materi¬ 
alize  in  this  instance.  The  problem  of  the  clamped  arch,  however, 
besides  being  of  some  practical  interest,  does  serve  as  a  very 
good  example  to  discuss  the  significance  of  the  energy  load  in 
relation  to  the  upper  buckling  load,  geometric  imperfections,  tran¬ 
sitional  buckling  modes,  and  varying  degree  of  approximation  of 
solutions.  The  clamped  arch  problem  also  contains  a  number  of 
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similarities  to  that  of  the  uniformly  loaded  spherical  cap.  The 
trend  of  results  of  the  arch  problem  serves  as  a  good  indicator 
of  the  behavior  of  the  latter  more  complicated  structure. 

Prom  the  experimental  side  it  is  possible  to  obtain  better 
control  with  fully  clamped  end  conditions.  This  is  especially 
true  at  high  temperatures  where  binding  at  the  supports  would  be 
a  problem, and  an  important  part  of  the  general  program  was  the 
Investigation  of  snap  buckling  due  to  high  temperature  creep. 

The  elastic  buckling  experiments,  which  were  performed  at  room 
temperature,  were  concerned  with  obtaining  complete  load  deflec¬ 
tion  curves,  including  the  unstable  regions,  for  a  wide  range 
of  arch  geometries, 

(a)  General  Formulation, 

The  circular  clamped  arch  loaded  with  a  single  central  con¬ 
centrated  load  as  shown  in  Pig,  10  is  analyzed.  With  the  usual 
assumption  for  thin  curved  beams,  the  axial  strain  e  is 

e  =  (1/R)  (ua-  w)  +  (1/2R2)  (wa)2  (19) 

and  x,the  change  in  curvature, is 

x  =  (1/R2)waa  (20) 

where  u  is  the  tangential  and  w  the  radial  displacement  function, 
R  is  the  radius  of  the  arch, and  a  is  the  polar  angle  measured 
from  the  center  of  the  arch.  Differentiation  with  respect  to 
a  is  denoted  by  the  subscript  a  , 

The  strain  energy  due  to  the  axial  deformations  ,  non- 
dimensionalized  by  division  by  the  factor  EtfR,is 


(21) 


Um  =  (1/2)  [e2da 

"P 


and  that  due  to  bending  Is 

ub  =  (t2M) 


P 

'  2. 
h  da 

■p 


(22) 


where  t  is  the  thickness  and  f  the  width  of  the  beam,  E  is 
Young's  modulus  and  p  half  the  Included  angle  of  the  arch.  The 
change  in  potential  energy  of  the  load,  in  non-dimensional  form, 
is  given  by 


Up  =  -  [P/EtRf  ]  (w)0_  0 


(23) 


where  P  is  the  concentrated  central  load.  The  total  energy  H  is 
given  by 


H  ■  um  +  “b  +  up  W 

The  tangential  displacement  function  u  can  be  eliminated 
from  the  energy  expression.  The  variation  of  the  total  energy  H 
with  respect  to  u  must  be  zero  for  equilibrium.  Since  u  only 
appears  in  Um,  the  membrane  strain  energy,  this  yields  immediately 

ea  =  0 

or 

e  =  C  (25) 

where  C  is  a  constant.  Substituting  in  (25)  from  the  strain  dis¬ 
placement  relation  (19)  and  rearranging 


Ua  =  w  -  (1/2R)  (Ma)2  +  CR 


Integrating  between  a  =  (3  and  a  =  -p  gives 

P  P  2  P 

[u]  =  J  [w-(l/2R)(wa)e;]da  +  [CRa] 

-P  "P  -P 

The  boundary  conditions  on  u  are 


u  *  0  at  a  =  +  p 


(26) 


The  left  hand  side  is  therefore  zero.  So’ v Inn;  th9  remainder  for 
the  constant  C  which,  from  (25)  is  equal  to  e,  yield3 


e 


lw-.(l/2R)(wa}'?]da 


-P 


(27) 


Substituting  tMs  value  for  into  the  energy  expression 
(24),  and  for  x  from  (20),  gives  the  total  energy  K  as  a  function 
of  w  only 


H  =  -V-'  |  ( w 

o! 


JL_ 

2R 


w2)da 


(28) 


2  ? 

where  the  brackets  in  (w0)  and  (n.rJ'  are  left  out.  In  the 
following  when  this  is  dene  it  will  always  mean  that  the  differ¬ 
entiation  is  performed  first, 

( b )  Approxlma to  Solu t i on 

A  first  numeric el  approximation  can  now  be  obtained  by 


assuming  a  function  for  w.  Let 
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where  Is  a  symmetric  and  w2  an  anti-symmetric  function  In  5  • 
5  Is  the  ratio  a/P  and  A  and  B  are  the  amplitude  of  the  two 
deflected  shapes.  Substituting  this  value  w  in  the  energy  ex¬ 
pression  gives 

H  *  (IAj^PHBAC^  +  pBC2  -  (1/2RP)  (A2C3  +  2ABC^  +  B^)]2 

+  (t2/24R^p3)  (A2C6  +  2ABCy  +  B2Cq)  -  (P/EtRf)  (AC9  +  BC^) 

(30) 


where  the  constants  are  given  by  the  integrals 


°i  ■  J  V5 

-1 


C2  =  w2dC 


C3 


2  .  r 


V  f 

-1 


c5 


°6  * 


1 

n  p 


(31) 


C7  = 


whsw25^ 


c8  = 


2  j  r 

w2« 


C9  =  [wx] 


5  =  0 


C10  =  ^W2^ 


5  =  0 


Introducing  non-dimensional  amplitudes 


a  =  (A/RP2) 


b  =  (B/Rp2) 


(32) 


2k 


and  load 

P*  =  (PR/Et2fp)  (33) 

and  a  geometric  parameter 

X  =  (p2R/t)  (3U-) 

the  energy  can  be  written  as 

H*  =  HX/p5  -  (X/4)  [a^  +  bC2  -(1/2)  (a2C3  +  2abC^  b2C^)]2 

+  (1/24X)  (a2C6  +  2abC?  +  b2CQ)  -  P*(aC9  +  bC1Q)  (35) 

Becai  se  of  the  symmetry  and  anti -symmetry  of  and  w2  It  ' 
can  be  seen  Immediately  that 


C2  =  C4  =  C7  =  °10  =  0 


The  energy  H  therefore  reduces  to 


H*  =  (X/4)[a01  -  (l/2)(a2C3  +  b2C5))2+  (l/24X)(a2C6  +  b2C0)-P*aC9 

(36) 

For  equilibrium: 


ft  # 

3H  _  n  3H  _  0 
■55  -  0  5F  "  0 


From  the  above  this  yields 


P*C9  =  (X/aHaCj-  (l/2)(a2C3  +  b2^))  (Cj-^a)  +  (l/l2X)(aC6)  (37) 


(bC^X/aHaC^  -  (l/2)(a2C3  +  b2C^)]  -  (l/l2X)bC8  =  0  (38) 


The  last  equation  has  the  solutions 


b  =  0 


(39) 
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b  =  (Z/C^ltiC^  -  (a2C3/2)  -  (C8/6X2C5)]  *  (40) 

If  the  quantity  under  the  root  is  negative,  the  only  real  solution 
is  b  =  0.  The  deformation  will  then  by  symmetrical*  The  condi¬ 
tion  that  a  real  solution  for  b  exists  other  than  b  *  0  is  there¬ 
fore  that 

aCx  -  (a2C3/2)  -  (Cq/Sx2^)  ±  0 

Which,  when  solved  for  a,  gives 

a  >  (Cj/C-j)  -  C  (C1/03>2  -  (Cq/IC^X2)]*  (14) 

a  <  (C^)  +  ((C^)2  -  {Cq/SC3C$X2))*  (42) 

For  "a"  real,  these  relations  will  be  satisfied  over  a  f lnite 
interval  if  the  quantity  in  the  root  is  greater  than  zero*  The 
necessary  condition,  therefore,  that  the  anti-symmetrio  component 
be  non  zero  is  that 

X2  2  c5  (43) 

If  condition  (43)  is  not  satisfied  the  load  deflection 
curve  is  obtained  by  setting  b  =  0  in  the  equilibrium  equation 
(37)  thereby  obtaining 

P*  =  (C6/l2C9\)a  +  (C2X/4C9)a[a-(C1/C3)](,a-(2C1/C3)]  (44) 

Such  a  curve  Is  plotted  in  Pig,  11,  It  c«i  be  considered  to  be 
composed  of  a  straight  line  due  to  the  bending  stress,  and  to  a 
curve  anti  symmetrical  about  a  =  Cx/C3  due  a**6!  stress. 
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This  is  exactly  the  same  kind  of  curve  as  was  obtained  for  the 
model.  The  energy  load  can  therefore  be  obtained  from  inspection 

as 

P*e  =  (1/12X)  (CjC^C^)  (1+5) 

The  range  of  stability  can  be  found  from  the  second  varia¬ 
tion  of  the  energy.  Substituting  b  =  0  into  the  total  energy  (36) 
and  differentiating  twice  with  respect  to  a  yields 

p  ^ 

=  (XyfeH^-aCj)2  -  aC^  +  (a2C2/2)]  +  (C6/l2  )  (1+6) 

For  stability,  this  second  variation  must  be  greater  or  equal  to 
zero.  Solving  for  a  gives 

a  <  (C^)  -  [(1/3)(C1/C3)2  -  (C6/9  C32X2)]i 
or 

a  1  (C1/°3)  +  IC1/3(C1/C3 )2  -  (C6/9  C32X2)]* 

These  portions  of  the  curve  are,  of  course,  the  portions 
with  a  positive  slope  in  Fig,  n.  The  region  of  the  load  deflec¬ 
tion  curve  that  is  unstable  vanishes  when  the  quantity  under  the 
root  in  (1+7)  and  (1+8)  becomes  zero.  An  unstable  region  will 
therefore  exist  for 

X2  >  C6/3C12  (1+9) 

and  instability  will  occur  after  the  equality  sign  in  (1+7)  is  satis¬ 
fied,  Substituting  this  value  for  a  into  the  load  deflection 
relation  (1+1+)  gives  the  upper  buckling  load  for  the  symmetrical 


(1+7) 

(1+S) 
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mode. 


C1C6 

T?CJC^  + 


i  °i 2 

3(^>- 


°6  ,3/2 


(50) 


If  the  Inequality  (43)  la  satisfied,  b  is  not  identically 
zero  but  is  given  by  the  expression  (4.0)  in  the  interval  given  by 
(4-1)  and  (42).  Substituting  for  b  from  (4-0)  into  (37)  yields  the 
load  deflection  relation  in  this  interval, 

P*  =  (C3C8/12\C£C9)[  (C^C-j)  -  a]  +  ( l/l2\ )  ( C6a/C9 )  (5D 

This  is  a  straight  line  passing  through  the  energy  load  P  s  P# 
and  a  *  C^/C^,  The  total  load  deflection  curve  will  be  then  be 
given  by  (44-)  except  for  the  values  of  a  in  the  interval  (41) 
and  (4-2)  where  it  will  be  given  by  (5D.  This  is  plotted  in  Pig. 11. 
The  value  of  a  at  which  the  non- symmetrical  transition  mode  comes 
in,  a^,,  is  obtained  from  the  equal  sign  in  the  inequality  (4-1). 
However,  buckling  may  already  have  taken  place  in  the  symmetrical 
mode  before  this  value  of  a  is  reached.  The  condition  for  buokling 
in  the  non- symmetrical  mode  must  be 

-  aU 

where  a^  is  the  deflection  at  the  upper  buckling  load  (Pig,  11), 
Substituting  for  a^  and  a^,  from  (4-7)  and  (41 ),  upon  making  both 
of  these  inequalities  into  equalities,  gives  for  a^,  <  a^, 

X2  ^  (52) 


The  behavior  of  the  arch  for  different  values  of  \  is 
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summarized  in  Fig,  12,  The  non- symmetrical  buckling  load  PT  la 
obtained  by  substituting  for  a  in  (44)  from  the  equality  condi¬ 


tion  of  (41) 


P#  = 

*T 


°1°6 


1  1  f°3°8 

12fl3«9  X  +  T3X[^ 


0-  0,  2 

A t iA 

°3 


8  * 


3030sV 


(53) 


The  stability  of  the  anti-symmetrical  mode  can  be  deter* 
mined  from  the  slope  of  the  load  deflection  relation  (5D*  This 
slope  is 

§5  =  (1/12\)[-(C3Cq/C5C9)  +  (C6/C9)] 
which  is  negative,  and  hence  unstable,  for 

C3C8  >C5°6  W 

The  first  term  in  (53)  is  the  energy  load  (45),  and  therefore 
from  (54)  and  (53)  the  anti-symmetrical  mode  will  be  unstable 
when 


This  result  concurs  with  the  general  conclusion  obtained  in  the 
previous  section, 

(c)  Determination  of  Buckling  Loads  by  Classical  Eigenvalue 
Theory. 

The  instability  of  the  arch  can  also  be  approached  by  con¬ 
sidering  small  displacements  about  an  equilibrium  position  as  in 
the  classical  eigenvalue  formulation.  The  buckling  load  is  deter¬ 
mined  in  this  procedure  as  the  load  for  which  geometrically  adja¬ 
cent  equilibrium  states  could  exist  in  addition  to  the  original 
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state  at  the  same  load.  Both  the  upper  buckling  load  Py  and  the 
load  f or  the  onset  of  transitional  buckling  can  be  determined 
by  this  method. 

Denoting  the  Initial  state  by  superscript  o  and  the  buckling 
displacements  and  strain  by  superscript  t  ,  the  change  In  strain 
energy  during  buckling  is 

Um  *  (V2)  1  *e'2  +  2e,e°)da 

-P 

P 

U£  =  (t2/2ij.)  f(x'2  +  2x°x»)da  (55) 

-P 

U»  *  -(P/EtRf)  [w»)o.  Q 

For  equilibrium  In  the  buckled  state,  the  variation  of  the 
total  energy  must  be  zero.  The  variation  with  respect  to  the 
tangential  buckling  displacement  u*  gives  simply 

e’  =0 

a 

or  by  integration 

e'  *  C«  (56) 

The  buckling  strain,  which  is  the  difference  between  the 
final  strain  and  the  strain  In  the  unbuckled  state.  Is  therefore 
a  constant.  From  the  definition  of  the  buckling  strain,  the  re¬ 
lation  between  that  strain  and  the  displacements  is 

e ‘  =  (l/R)(u»  -  w')  +  (l/2R2)(w»2  +  2w°w»  )  =  C« 
u  a  a  a 


(57) 


30 


whioh,  upon  solving  for  u*  gives 

ct 

u*  *  w'  -  (l/2R)(w»2  +  2w°w»)  +  C»R  (58) 

01  0 1  u 

Integrating  between  a  *  (3  and  o  s  -P  and  utilizing  the  clamped 
boundary  conditions  (u  =  0)  yields 

P  P 

C*  *  -(1/2RP)  jCw*  -(1/2R)  (w*^  +  2w®w»I  )]do  (59) 

"P 


Substituting  this  value  for  C*  into  (56),  the  change  in  energy  due 
to  buckling  can  be  written  as 


H»  «  (iAr2P)J 

I 

-3 


P  p  > 

p Tw*  -(l/2R)(w^  +  2w®w^  ))dO  , 

"P  J 

Pf  P 

o. 


+  (t2/2UR^)  J(w«a  )2do  -  Jie°/2RP)  J[w»-(l/2R)  (w^+2w°w»  ))do^da 


-P 


P  O  ^ 

+  (tV12R  )  H°w*  da  -  ( P/EtRf )  ( w  * ) 
i  aa  < 

-P 


a*  0 


(60) 


Because  the  initial  state  is  an  equilibrium  state  all  the 
linear  terms  in  w  must  cancel  out  either  algebraically  or  due  to 
boundary  conditions.  If  it  is  assumed  in  addition  that  w  is  suffi¬ 
ciently  small  so  that  terms  of  higher  order  than  second  can  be 
neglected  (in  accordance  with  the  basic  assumption  of  classical 
buckling  theory),  then  the  energy  change  is  reduced  to 
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H»  =  (l/4R2p) 


p 


1 


p 


[w*  -  (l/R)w°w*a)]da',»  +  (t2/2lj.R^)  fw*.2  da 


-P 

+  (l/2R2)e° 


w*2da 

a 


aa 


(61) 


The  pre-buckling  strain  e  can  be  taken  outside  the  integral  since 
it  is  a  constant.  Assume  that 

w1  «  A*w^  +  B’Wg 

where  is  a  symmetric  furotion  in  a  and  v2  is  ant  1- symmetric. 
The  functions  and  w2  are  therefore  orthogonal  in  the  interval 
-  P  <  a  <  P  .  Substituting  into  the  energy  expression  (61)  gives 


H*  =  (l/l4irp)<|  |  [ A,w1  ♦  BW,  -  il/RH0(A»wla  +  B'w2a))da 


+  (t2/24R]+)  f  (A'wlaa+  3»w2aal2da  +-(l/2R2)s°  J(A'wla+  B'w2a)2da 

-P  -P 


(62) 


Because  of  the  orthogonality  properties,  this  reduces  to 


H 


2' 

«  =  A1  > 


(1/1^2P)[  f(w1  -  (1/R)w°wla)da]; 


"P 


+  (t2/2l^)  jWl2ada  +  ( e°/2R2)  f wn2da 


r 


+  B»2',  (t2/2l4^^) 


-P 

0 


U  1 

la 

4 

P 


w2aada  +  (e°/2R2)  |'w2a  da 


-P 


(63) 
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For  equilibrium  In  the  buckled  state 


dH' 

"35* 


0 


Each  of  the  two  expressions  In  parentheses  In  (63)  must  therefore 
be  zero. 


From  the  second. 


This  gives  the  strain  In  the  arch  when  the  non-symmetrie 
mode  buckling  load  Is  reached.  From  the  form  of  the  expression 
(65)  It  Is  seen  that  the  stress  Is  equal  to  the  second  mode  Euler 
load  of  a  clamped  column  of  the  same  length  as  the  arch.  By 
substituting  from  (27)  for  e°  the  deflection  at  buckling  can  be 
obtained 

P  A 

(1/2RP)  |  (v°-(l/2R>H°2]da  =  (t2/12R2p2)  w2^dC/ (66) 

-P  -1  /  -1 

Assuming  that  w°  =  Aw^(£)  and  introducing  the  non-dimen¬ 
sional  amplitude  a  defined  in  (32),  and  the  parameter  \  defined 
by  (3(4.)  the  above  expression  becomes 
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wx25d5 


-  2a  +  (1/3X2) 

-1 


-1 


1  2  ,r 
w2^d^ 


=  0 


Upon  substituting  for  the  integrals  from  (31), 
a2  -  (2C1a/C3)  +  (Cq/SX2^^)  =  0 

This  is  the  same  condition  as  obtained  previously,  (l+l),  for  the 
amplitude  a^,  at  which  the  transitional  mode  appears.  . 

The  first  equilibrium  equation,  obtained  by  setting  the 
variation  of  H',  (63),  with  respect  to  A' equal  to  zero,  can  be 
written,  using  (27)  for  eQ, 


(1/20) 


f 

vj 


(w1  -  (l/R)w°wla)dcc 


+  (t2/l2R2) 


wlaada 


-  (1/20R) 


f(w°  -  (l/2R)w°2)da 

-0 


s 

,  1 

“la  da  j 

_ 

W 

"  J 

Again  assuming  that  w°  =  Aw^  the  equation  becomes 

(1/2  3) 


r 

p 

2  0  D  P 

■  1 

■'[w1  -  (A/R)w1^]da 

+  (t2/l2R2) 

J 

P 

-B 

w,  __da 
laa 


5 

fi 

P 

j 

'[Awx  -  ( A2/2R)w^2) da > 

J 

'wladaf 

*  1 

P 

L-1 

-  (1/2R0) 


Using  the  previous  notation,  this  equation  reduces  to 
a2  -  (2G1a/03)  +  (2/3  )  ( Cj/C^  )2  +  (C6/9X2C32)  =  0 


(67) 


(68) 
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which  is  the  same  condition  as  obtained  previously  for  the  deflec¬ 
tion  a^j  at  the  upper  buckling  load  (47  )•  The  strain  at  which  this 
occurs  can  then  be  determined  from  (27)*  The  corresponding  load 
Pjj  can  be  obtained  by  setting  a  =  a^  in  the  load  deflection  rela¬ 
tion  which  would  be  identical  to  (44)* 

This  analysis  shows  that  although  the  problem  is  basically 
a  nonlinear  one,  both  the  upper  buckling  load  for  symmetrical 
buckling  and  the  non- symmetrical  transitional  buckling  load  can 
be  determined  by  the  classical  eigenvalue  procedure.  The  essential 
difference  between  this  problem  and  the  standard  ones  is  that  non¬ 
linear  strain  displacement  relations  must  be  used  in  describing 
the  pre-buckllng  states, 

(d)  Analysis  of  a  Geometrically  Imperfect  Clamped  Arch. 

It  is  of  Interest  to  compare  the  behavior  of  the  clamped 
arch  with  initial  imperfections  with  the  observations  in  the 
preceding  general  discussion.  It  will  be  assumed  that  the  devia¬ 
tions  from  the  perfect  circular  shape  are  given  by  the  functions 
we  and  u®,  The  elastic  membrane  strain  can  then  be  written  as 

e  =  (l/R)(Ua  -  «)  ♦  (l/2R2)(w£  +  2waw®)  (69) 


and  the  change  in  curvature,  as  before,  is  given  by  (20).  The 
energy  expressions  (21)  (22)  and  (23)  still  hold.  By  varying  the 
total  energy  with  respect  to  the  displacement  u,  as  previously, 
the  axial  strain  is  found  to  have  a  constant  value  given  by 


-  (1/2RP) 


f[w  -  (l/2R)(w2  +  2waw®)]da 


4 


e  = 


(70) 
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and  therefore  the  total  energy  can  be  expressed  as 


H  =  (lAR2p) 


f[w-(l/2R)(w2  + 


2waw®)]da 


!>  +  (t2/2^R2) 


4 


2 

w _ da 

aa 


-  (P/EtRf)  (w)a=  Q  (71) 

It  Is  assumed  that  w  Is  given  by  (29)  and  the  initial  imperfections 
by  we  =  Dw^(0  +  Sw2(0  where  w^  and  Wg  are  the  same  functions 
as  in  (29).  The  total  energy  can  then  be  written  as 

H*  *  (XAljaC^  +  bCg  -  (l/2)(a2C3  +  2abC^  +  b2C^) 

I2 

-  [adC^  +  (bd  +  ae)C^  +  boC^]| 

+  (l/2i4.X)(a2C6  +  2abC_  +  b2CQ)  -  P*(aC9  +  bC10)  (72) 


where  the  quantities  a,b,d,e  are  non-dimensional  amplitudes  obtained 
by  dividing  the  corresponding  capitalised  quantities . by  the  factor 
Rp2,e.g.  (32), 

Because  of  the  orthogonality  of  the  functions  w^  and  Wg, 

c2  =  =  c7  =  C10  =  0 

H*  therefore  can  be  reduced  to 


H*  =  ^/^[a^-d/sKa2^  +  b2C^)  -  (adC3  +  beC^)]2 

+  (lM\)(a2C6  +  b2CQ)  -  P'*aC9  (73) 

For  equilibrium,  ^  =  0,  =  0, 
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which  results  in  the  following  two  equations: 

-  (VsHaC^-d/aMa2^  +  b20^)-(adC3  +  beC^)] (Cj-aCyd^ ) 

+  (l/l2\)(aC6)  (74) 

(X/^HaCj-dyteJa2^  +  b2C^)-(adC3  +  beC^JHbC^  +  eC^] 

-  (l/l2X)bC8  =  0  (75) 

Unlike  the  corresponding  equation  for  the  geometrically 
perfect  system,  (38),  b  =  0  is  no  longer  a  solution.  The  deform¬ 
ations  in  the  anti-symmetrical  mode  therefore  commences  immediately 
at  the  onset  of  loading. 

Upon  changing  the  origin  in  the  load  deflection  plane  by 
the  transformation 
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Prom  (78)  it  la  aeen  that  the  amplitude  of  the  anti¬ 
symmetric  mode  b  is  independent  of  the  sign  of  a’,  i.e.  b  is 
symmetric  about  a»  *  0.  The  load  deflection  relation  (77)  consists 
of  a  membrane  contribution  (the  quantity  containing  \ )  and  a  bend¬ 
ing  part  (containing  l/\).  The  bending  term  is  linear  in  a*  where- 

a  which 

as  the  membrane  term  ls/oubic  in  a*  /  is  anti -symmetrical  about  and 
passes  through  a'  =  0.  The  load  deflection  curve  is  therefore  com¬ 
posed  of  a  straight  line,  and  a  curve  anti-symmetrical  about 
a  =  [(C^/C^-d],  The  energy  load  is  therefore  given  by 

Pe  *  (C6/12xC9>[(Gl/C3>-d]  (79) 

If  the  initial  imperfections  are  anti-symmetrical,  d  =  0,  and  the 
energy  load  is  the  same  as  for  the  geometrically  perfect  arch  (45). 

The  behavior  of  the  clamped  arch  is  therefore  similar  to 
that  of  the  model  in  that  the  behavior  in  the  presence  of  anti- 
symmetrical  imperfeotlona  is  identical  to  the  type  A  imperfection 
of  the  model.  Symmetrical  Imperfections  of  the  arch  would  be  of 
type  B  and  correspond  to  a  change  in  the  overall  geometry  and 
would  result  in  a  change  (lowering)  of  the  energy  load  as  noted 
in  (79). 

(e)  More  Exact  Solution  of  the  Nonlinear  Deformations  of  the 
Clamped  Aroh. 

A  better  approximation  to  the  load  deflection  relations  and 
to  the  buckling  loads  of  the  arch  can  be  obtained  by  including  more 
terms  in  the  expression  for  the  deflected  shape.  Symmetrical 
deformations  will  first  be  considered.  The  deflected  shape  is 
assumed  to  be  given  by 
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W=S=1  .Vn(S> 


(80) 


Ajj  Is  then  the  amplitude  of  the  deflection  function  wn  and  K  =  ®/P 
as  before.  Substituting  the  deflected  shape  (80)  into  the  expires* 
sion  (28)  for  the  total  energy  yields 


fl 


H*  =  (\A)j 


[  |  nnwn)d?  -  * 


1 

i1 


r 

Z 


Vat1 


+  (l/2llX) 


[  Z  a„w, 


]2<U  - 


P*[ 


r 

2  awl 
1  n  n  £  = 


(81) 


where 


a 


n 


and 


If  wnf.  and  wn^  form  orthogonal  sets  such  that 


(82) 


(83) 
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The  total  energy  Is  reduced  to 


/• 

r 

r  ,  1 

■ 

*  0 

=  (Vfc>- 

2  [a 

L1  n-i 

V?1  -  *  ? 

1  -1 

”n5d51- 

»  £  ,  2, 


+  U/24X)  I  [<.  J  w^dtf  -  P*  Z  i^n\m  0) 


(84) 


The  equilibrium  condition  la 

c*° 

which  leads  to  r  equations  of  the  form 

(X/2) 


m 

1 

r 

r  o  1 

2  [an 

1  n  J 

k.. 

wndQ  "  *2  [an  J 
L  1  -] 

> —  - 

X 

3 

/n 

ffc 

' - V 

/ 

■> 

L  1 

■ 

\ 

|Hnd<  -  Si  ,[ 

>  ♦  (l/l2\)a 

11  V 

c1 

L  -1 

s< 

"A 

fXl  *  0  (85) 

n5  =  0 


The  following  notation  is  Introduced: 


'In 


-1 


»nd5 


°2n  ' 


(86) 


'3n 


V  "  twn\=  0 


The  equilibrium  equations  for  symmetrical  deformation  of  the  arch 
can  then  be  written  as 
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<x/2)l 


I  Si0!!!"  *  | 


4°2n 


1 

°ln-  Van  (87) 


J 


The  criterion  obtained  previously  for  the  occurrence  of 
antl-symmetrlcal  buckling  (66)  remains  essentially  unchanged.  The 
only  modification  would  be  that  the  more  general  expression  (80) 
would  be  used  to  represent  the  pre-buckled  state.  This  generali¬ 
zation  of  (66)  leads  to 


2  [a, 


n 


vndO  -  (1/2)  2  [a2 


w2  ] 
wn£J 


K/6X2 


(88) 


where 


(89) 


and  wa  is  the  antl-symmetrlcal  buckled  shape.  Upon  using  the 
expressions  (86)  for  the  constants,  the  criterion  for  antl- 
symmetrlcal  buckling  becomes 

f  V=ln  -  ftnC2n  *  (90) 

To  obtain  the  antl-symmetrlcal  mode  buckling  load  the  set 
of  equations  (87)  has  to  be  solved  simultaneously  with  (90), 
Substituting  (90)  into  (87)  yields 

(K/l2X)(0ln  -  an02n]  +  U/l2ManCJn  -  P*C^  =  0  (91) 

which  can  be  solved  for  a_ 


an  88  <12XP\n  ‘  Kcln)/^C3n“  Kc2n> 


kl 


(92) 

This  value  for  an  can  now  be  substituted  back  into  the  buckling 
criterion  (90 )# 

(12XP*)2  |  [C2n<^n/(c3n  -  K02n)2]  -  <24xp£>  il  c1Iciaclv/{  c3n.KC2n,2, 

+  (K/3X2)  +  K  2  [C^n(2C3n  -  KC2n)/(C3n  -  KC2n)2]  =  0  (93) 

This  can  be  written  as 

(12  XP*)2^  -  2!|AP*D2  +  (K/3X2)  +  KD3  *  0 
and  therefore  the  asymmetrical  mode  buckling  load  Is  given  by 
Pj  *  ( l/l2  X)  <  (Dg/I^)  +  [  (D2/D^)2  -  (  KI>3/d^  )  -  (K/3X2D1)]^>  (914-) 

where 

°1  =  |  t°2n°4i/(03n  *  K02n)2> 

B2'f  t°lnc3nV(03n  -  KC2n)2’  <’5) 

°3  ‘  f  "4<2C3n  *  KC2n>/(C3n  -  K02n>2> 

The  seoond  buckling  mode  can  therefore  be  determined  fairly 
easily  by  evaluating  the  C’s  and  D»s  and  solving  a  quadratic  for 
12\P*. 

For  a  clamped  arch  a  general  displacement  function  that 


satisfies  the  boundary  conditions  Is 
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wn  *  id  +  cos  n  *53#  a  odd 

(96) 

wn  *  i(l  -  oob  n  icOi  n  even 

This  will  also  satisfy  the  orthogonality  relations  (63).  The 
constants  (86)  oan  then  be  readily  evaluated. 


°2n  * 

C3n  * 

m  lt  n  odd 

0#  n  even 

The  quantity  K  Is  the  coefficient  In  the  equation  for  the 
second  mode  buckling  load  for  a  damped  column,  namely, 

K  ■  20,16  or  2.0426it2 

The  constants  (95)  then  beoome 

*  2  (4/n2n^(n2-2#0426)2],  n  odd 

*  0,  n  even 

D?  *  S  I4/^(n2-2.0426)2],  n  odd 
d  1 

“0,  n  even 

D,  *  E  [4(2n2-2. 0426  )/n2^(n2-2, 0426 )2],  all  n 
5  1 


(98) 


The  values  of  these  constants  are  for 


r  *  Is 

Dx  »  4  x  0.9210/rc6,  D2  SU  0# 9210/i^,  Dj  *  -  4  *  0.00387A14- 
r  =  2: 

and  Dg  unchanged  and  *  4  x  0.3498/x^ 

For  r  *  6,  which  is  sufficient  to  determine  the  constants 
to  3  figures  after  the  decimal  point,  the  results  are 

Dx  =  4  x  0.9233 A6,  D2  -  4  x  0.94.36/it^,  D3  =  4  *  0.3960/^ 

The  antl-symmstrlcal  transitional  buckling  load  can  then 
be  readily  obtained  from  (94)  for  these  three  cases 


r  =  1 

P*  =  (x2/l2X) 


f 

1  +  [1.0858  -  (0.l847*Vx2]* 


L. 


r  =  1 


r  = 


(x2/l2X) 
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1  +  [0.2245  -  (0.l847iA/X2)]4 


K  -  ( 


tc2/12\)]  1.0219  +  [0.1685  -  0.1843^/X2]^ 


.1 


l 
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These  three  functions  are  plotted  in  Figure  13.  There  is 
an  appreciable  difference  between  the  curves  for  1  and  2  terms, 
but  not  much  between  the  curves  for  2  and  6  terms. 

Determining  the  upper  buckling  load  for  the  symmetrical 
mode,  Py,  is  somewhat  more  involved.  It  is  therefore  only  attempted 
for  r  =  1  and  r  =  2,  For  r  =  1  the  upper  buckling  load  is  obtained 
immediately  by  substituting  the  constants  (97)  into  the  expression 
for  one  term  obtained  previously  (50)  observing  that 


Cln  *  °l»  °2n  s  c3»  °3n  =  °6»  °4n  *  c9  when  n  e  1 


This  gives 

3/2 

pj  =  (*?/l2X)  +  (Xi^A)  C (U./3W+)  -  (1/9 x2)]  (99) 

The  energy  load  for  a  single  term  is  obtained  by  substitut¬ 
ing  the  values  of  the  constants  into  (45)  thereby  obtaining 

P*  *  *?/l2X 

The  next  approximation,  r  **  2  will  be  considered, 
equilibrium  equations  187)  it  is  seen  that 

p*<w<vVi2>->  _ 

"In  ■  n  2n  Clm  *  am°2m 

which  can  ba  solved  for  a  In  terma  of  a^  and  the  load  P* 

.  P*<Chnel.-clnlU-»mChnCa.)  *  «m>clnCl./12X> 

^  *P*°2n%»  +  <»n/12W(c2n03n-C3n02n>  +  <  Clmc3r/12X) 

For  r  =  2  this  reduoea  to 

x2a2  =  -  (kx  -  P*)/(k2  +  P*) 

where 

\  =  a/48jj.,  k2  =  l/p.[  (a/l6 )  -  (1/3)1, 

2  2 
a  =n  a^,  and  jt  =  X/x 


(101) 


(102) 


(100) 
Prom  the 


The  first  of  the  equilibrium  equations  (87)  can  then  be  written 


as 


[Ax  +  “  (a^/8)  -  a|«k/2)l  [1  -  (a^A)] 


+  (a17:D/2l4.X.<:)-(2n2pVx)  *  0 


(103) 


Upon  substitution  for  n^ag  from  (101)  and  the  use  of  the  notation 
(102),  a  relation  between  a  and  P*  can  be  obtained.  This  is 


( a~(a/8) ]  (kg  +  pV-O^-^X^  +  P*)-(l/2)(lc1-P#) 
+  ( ( a/24. ^ )  -  ( 2P*/ |i)  ]  ( k~  +  P*)2  »  0 


[l-(aA)l 


(104) 


The  load  deflection  relation  for  symmetrical  deformation  is 
given  explicitly  by  this  expression.  Because  of  the  complexity  of 
the  expression,  the  load  deflection  curve  can  be  best  determined 
by  first  introducing  new  variables  a  and  p  defined  as  follows: 

a  =  a  -  4  p  =  4®^?  ”  4  (105) 

Equation  (104)  can  then  be  written  as 

Ko  +  Kxp  +  KgP2  -  4p3  =  0  (106) 

where 

Kq  =  a3 [36  +  3u2(9a2-  116)] 

Kx  =  a 2 [ - 1 2  +  3p2(6a2  -  120)]  (107) 

K2  =  a[-20  +  3 ix2 (a2  -  20)] 

The  solution  to  (106)  is  anti symmetrical  in  p  about  the 
origin  a  =  p  =  0,  The  energy  load  is  therefore  Immediately 


obtained  as 

Pe  a  0 

or  from  (105)  and  (102) 

P*  «  1^/12 

The  energy  load  is  therefore  unchanged  by  going  from  1  to  2  terms 
in  the  expansion  for  the  deflected  shape.  The  addition  of  the 
second  term  is  found  not  to  alter  the  relation  between  the  load 
and  the  bending  deformation  and  the  symmetry  of  the  membrane 
stresses.  The  influence  of  the  seoond  term  is  therefore  similar 
to  that  of  a  type  A  imperfection.  This  result  would  probably  not 
be  true  by  the  inclusion  of  more  terms  in  the  deflected  shape 
function,  but  it  is  believed  that  the  consideration  of  additional 
terms  will  not  have  any  significant  effect  on  the  energy  load. 

The  upper  buckling  load  for  the  symmetrical  mode  is  most 
easily  obtained  by  plotting  p  against  a,  and  determining  the  valve 
of  p  at  the  limit  point  from  the  resulting  plot.  This  value  can 
be  transformed  back  to  the  original  coordinates  using  (105)  and 
(102),  Performing  this  procedure  for  a  range  of  the  geometric 
parameter  X,  *  plot  of  P^  against  X  is  obtained.  This  is  shown 
in  Figure  13  along  with  the  energy  load,  the  one  term  approxima¬ 
tion  for  Py,  and  the  transitional  buckling  load  P^  . 
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(f )  Experiments  on  Elastic  Snap  Buckling  of  damped  Arches. 

The  experimental  part  of  the  study  was  concerned  with 
obtaining  complete  load  deflection  curves  for  a  wide  range  of 
geometries  of  clamped  arches.  It  was  of  special  interest  to  obtain 
the  unstable  regions  of  these  curves  and  to  investigate  the  mode 
shapes  at  various  stages  of  loading.  The  experimental  setup  was 
constructed  so  that  creep  snap  buckling  tests  as  well  as  elastic 
tests  were  possible.  In  this  section  only  the  elastic  tests  con¬ 
ducted  at  room  temperature  are  considered.  Figure  14  shows  a 
schematic  drawing  of  the  apparatus  and  Fig.  15  and  16  are  photo¬ 
graphs  of  the  actual  apparatus. 

The  arch  was  clamped  between  two  fixing  blocks  at  each  end. 

It  was  necessary  to  make  a  series  of  these  blocks  of  different 
angles  to  cover  the  range  of  geometry.  The  fixing  blocks  were 
bolted  to  a  very  stiff  auxiliary  frame  which  rested  freely  upon  the 
main  support  frame.  This  was  necessary  to  allow  the  auxiliary 
frame  to  expand  freely  for  the  high  temperature  tests.  To  ensure 
complete  rigidity  of  the  clamping,  the  fixing  blocks  were  backed 
up  against  two  other  blocks  which  were  bolted  and  doweled  to  the 
auxiliary  frame. 

The  deformations  of  the  arch  were  measured  at  several  points 
with  dial  gauges.  These  were  arranged  outside  the  furnace  with 
extension  rods  resting  on  the  arch.  The  extension  rods  were  counter- 
weighted  with  springs  to  minimize  the  effect  of  their  weight  on  the 
behavior  of  the  arch.  Experiments  with  only  the  center  gauge  and 
all  the  gauges  confirmed  that  this  was  achieved. 


The  load  was  applied  at  the  center  through  a  knife  edge 
fixture  and  a  loading  rod  extending  out  of  the  furnaoe,  To  obtain 
the  complete  load  deflection  curve,  a  system  of  a  jack  and  a  strain 
gauge  instrumented  load  cell  in  series  with  the  applied  load  was 
utilized.  The  load,  which  was  greater  than  P^,  rested  upon  the 
jack  and  the  entire  load  deflection  relation  was  obtained  by  lower¬ 
ing  the  jaok  in  small  Increments  -  in  effect  a  deflection  controlled 
procedure  -  and  noting  the  load  supported  by  the  arch  from  the  load 
cell  and  the  corresponding  deflections  from  the  dial  gauges. 

The  specimens  were  made  from  2021J.T1J.  aluminum  alloy  and  were 
one  inch  wide,  three  sixteenths  Inches  thick,  and  nominally 
3I4.  Inches  between  supports.  They  were  rolled  into  a  circular  arc 
using  a  three  roller  sheet  metal  roll.  The  geometry  was  determined 
by  measuring  the  central  height  after  the  specimens  were  bolted  into 
the  auxiliary  frame;  The  clamping  resulted  in  a  slight  pre-stress- 
ing  for  some  of  the  specimens,  particularly  at  the  lower  values  of 
X  ,  as  it  was  very  difficult  to  get  the  angle  of  the  fixing  blocks 
to  coincide  exactly  with  the  base  angle  of  the  specimens. 

In  all,  1^  specimens  were  tested  for  X  ranging  from  3.69 
to  16,25,  Table  I  gives  the  geometries  and  the  buckling  loads  ob¬ 
tained,  All  the  experiments  were  completely  elastic  as  the  speci¬ 
mens  were  found  to  return  to  their  original  shape  to  within  a  few 
thousandths  of  an  inch  upon  unloading.  The  lower  buckling  load  of 
the  clamped  arch  is  greater  than  zero  and  therefore  the  arch  cannot 
remain  in  the  buckled  state  when  the  load  is  removed. 

The  experimental  buckling  loads  (the  maximum  Doints  in  the 
load  deflection  curves)  are  compared  with  the  theoretical  results 


in  Pig,  13#  Except  for  the  two  points  with  the  lowest  values  of 
X,  the  agreement  with  the  two  term  solution  is  very  good.  It  is 
probable  that  the  disagreement  in  those  two  cases  was  due  to  the 
clamping  stresses. 

Pig,  17  shows  complete  load  -  central  deflection  curves  for 
three  values  of  X  along  with  corresponding  theoretical  predictions. 
The  complete  theoretical  load  deflection  curves  (for  1  and  2  terms) 
are  shown  for  the  X  =  11,62  case.  Although  the  general  trends  and 
especially  the  buckling  load  values  are  in  close  agreement,  the 
curves  by  themselves  are  not  conclusive  as  to  the  buckling  mode 
in  the  unstable  region.  It  is  noted,  however,  that  for  large  X, 

Pg  is  appreciably  greater  than  P^,  and  the  experiments  agree,  as 
expected,  with  PT, 

Although  the  antl-syranetrlcal  mode  governs  the  buckling 
process  for  X  >  10.6,  the  actual  amplitude  of  the  anti -symmetrical 
mode  is  relatively  small  and  could  be  readily  overlooked  in  a  super¬ 
ficial  test.  The  shapes  at  buckling  for  various  arches  are  shown 
in  Pig,  l8.  The  Increased  waviness  in  the  shapes  for  the  larger 
X  values  and  the  corresponding  need  for  more  than  one  term  in  the 
theoretical  solution  can  be  readily  noted.  The  asymmetry  of  the 
shapes  can  be  noticed  for  the  large  X  values  but  is  very  small 
compared  to  the  overall  symmetrical  deformation.  Pig.  19  shows  the 
deflected  shapes  for  two  values  of  X  at  the  upper,  lower,  and 
middle  buckling  points.  Even  at  the  middle  buckling  position,  when 
the  anti-symmetrical  amplitude  is  largest,  the  asymmetry  is  small 
for  the  X  =  16,25  case  and  almost  non-existent,  as  expected,  for 
X  =  5.35.  These  results  show  that  although  the  buckling  of  a 


50 


system  may  be  governed  by  anti-symmetrlcal  modes,  these  modes  may 
not  be  readily  detectable  except  by  very  careful  experimentation. 
These  observations  are  believed  to  be  pertinent  to  the  related  prob¬ 
lem  of  the  buckling  of  spherical  caps, 

(g )  General  Discussion  of  Results  for  the  Clamped  Arch. 

1.  The  experimental  results  were  in  good  agreement  with  the 
general  theoretical  predictions  and  showed  that  for  large  X  buckling 
is  governed  by  anti-symmetrical  modes  -  the  buckling  load  corres¬ 
ponded  to  the  transitional  buckling  values  P^,.  The  amplitudes  of 
the  anti-symmetrical  modes  were,  however,  relatively  small.  As 
expected,  the  energy  load  was  a  lower  bound  on  the  experimental 
buckling  loads,  and  could  have  served  to  give  reasonable  predictions. 

2.  Although  a  change  in  buckling  mode  occurs  (theoretically) 
at  X  *  10.6,  the  transition  Is  a  smooth  one  and  the  buckling  load- 
geometry  (X),  curve  is  also  smooth.  The  behavior  is  therefore  un¬ 
like  the  sharp  divisions  obtained  for  plate  buckling,  and  is  a 
further  reason  for  the  difficulty  in  detecting  bucklir.g  mode  changes. 

3.  The  general  features  of  the  nonlinear  behavior  of  the 
clamped  arch  were  found  to  be  similar  to  those  of  the  simple  mech¬ 
anical  model  and  to  conform  with  the  general  results  obtained  for 
such  nonlinear  systems, 

Ij..  The  one  term  solution  gave  the  proper  overall  behavior  but 

was  appreciably  in  error  on  the  actual  magnitudes  especially  for 
the  higher  values  of  X  .  The  two  term  solution  gave  fairly  good 
results  -  the  additional  corrections  to  the  six  term  solution  for 
the  anti-symmetrlcal  buckling  load  are  relatively  small. 
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5*  The  energy  load  waa  the  same  for  both  the  one  and  two  term 

solutions.  Although  the  energy  load  would  not  be  unchanged  in  the 
complete  solution,  there  are  strong  Indications  that  the  energy  load 
Is  a  relatively  insensitive  function  of  the  buckled  shape. 
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Creep  Buckling  of  Arches. 

The  creep  buckling  problem  was  studied  as  an  Important  prob- 
lem  In  its  own  right  and  as  an  example  in  a  difficult  field  where 
the  energy  load  concepts  may  have  application.  The  possibility  that 
the  energy  load,  as  defined  for  the  initial  geometry,  may  be  a  lower 
bound  on  the  load  required  for  snap  buckling  in  the  presence  of 
creep  was  the  point  of  particular  interest.  However,  it  became 
evident  early  in  the  work  that  the  energy  load  based  on  the  initial 
geometry  was  not  a  lower  bound  in  the  above  sense.  The  only  signifi¬ 
cance  of  the  energy  criterion  in  creep  problems,  at  least  for  the 
model  and  the  aroh,is  that  the  "membrane"  forces  corresponding  to 
the  original  energy  load  must  be  exceeded  before  buckling  can  occur. 

Considerable  insight  into  the  snap  creep  buckling  problem 
can  be  obtained  from  the  analysis  of  the  simple  model  of  Fig,  2. 

A  related  study  was  performed  by  Hult  [18],  but  his  model  did  not 
contain  the  central  spring  which  would  correspond  to  the  bending 
rigidity  in  the  case  of  the  arch.  The  energy  load  for  Hult’s  model 
is  therefore  identically  zero  for  all  geometries  and  cannot  be  used 
to  investigate  the  significance  of  the  energy  load  in  creep  snap 
buckling  problems.  In  addition  Hult*s  model  is  statically  deter-  • 
mlnate  whereas  the  foroe  distribution  in  the  more  general  model 
depends  upon  the  relative  stiffnesses  and  the  creep  law  as  well  as 
upon  the  geometry. 

The  problem  of  oreep  snap  buckling  can  be  formulated  as  the 
determination  of  the  load  deflection  curve  for  superimposed  loading 
at  successive  stages  of  creep  of  the  system.  When  the  geometry  and 
force  distribution  is  such  that  the  load  increment  required  to  reach 
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either  or  P^  (whichever  is  lower)  becomes  zero,  then  snap 
buckling  will  ocour.  The  analysis  is  based  upon  assuming  that  the 
superimposed  load  defleotlon  curve  can  be  calculated  from  elastic 
theory.  The  incremental  stresses  and  strains  due  to  buckling  are 
taken  to  be  elastic.  This  assumption  agrees  with  the  commonly  used 
creep  laws  since  the  snap  buckling  process  is  relatively  rapid. 

The  equilibrium  equation  for  the  model  of  Pig,  2,  as  given 
by  (2),  can  be  written  with  slightly  greater  generality  as 

P  *  P  ♦  2T(a-6)/L  (108) 

where  P  Is  the  force  In  the  compressed  spring,  and  the 

other  symbols  are  as  defined  previously.  The  strain  In  the  arms 

is  given  by  (3).  In  the  following  analysis  It  Is  convenient  to 

oonslder  compressive  strains  as  positive.  Both  the  arms  and  the 

spring  are  considered  to  creep  under  the  applied  forces. 

The  model  Is  subjected  to  a  oonstant  applied  load  P  .  The 

o 

forces  In  the  arms  and  in  the  spring,  and  the  apex  deflection  at 
any  later  time  are  designated  TQ,  PQ  and  6Q  respectively.  These 
quantities  are  also  related  through  (108),  The  load  deflection 
relation  for  an  additional  load  P  can  now  be  obtained. 

The  total  force  In  the  spring  due  to  PQ  +  AP  is 

P  =  P0  +  k(6-60) 

where  6  Is  the  total  displacement  from  the  unloaded  position.  Prom 
the  strain  expression  (3),  the  total  force  in  the  arms  is 

T  =  TQ  +  (AE/2L2)(6-60)(2a-6-5o) 
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The  total  forces  must  also  satisfy  the  overall  equilibrium  equation 
(108),  Upon  substitution  of  the  above  expressions  for  P  and  T  into 
(108)  and  subtracting  the  equilibrium  conditions  for  the  original 
loading  state,  the  load  deflection  relation  for  the  additional  load 
AP  is  obtained  as 

AP  *  [k-(2T0/L)l(6-6e)  +  ( AE/L3 ) (a-b) (5-5Q) (2a-6-6o )  (109) 

AP  is  plotted  against  6  in  Pig.  20.  This  is  the  instanta¬ 
neous  load  deflection  curve  for  the  system  at  any  stage  of  creep. 
The  curve  again  consists  of  a  straight  line  and  a  curve  anti- 
symmetrical  about  6  *  a.  The  energy  load  is  therefore 

APe  =  [k-(2T0/L)](a-60)  (110) 

For  the  fully  elastic  case  this  result  agrees,  as  expected,  with 
that  obtained  by  subtracting  (2)  from  the  elastic  energy  load, 

P  =  ka,  (9).  The  energy  load  for  additional  loading,  AP  ,  is 
aero  when 


TQ  =  kL/2  (111) 

which,  from  (9)  and  (2),  is  also  the  value  for  the  force  in  the 
arms  (the  "membrane”  force)  at  the  elastic  energy  load  P  .  The 
condition  6  =  a  is  a  special  case  for  the  vanishing  of  APa. 

The  upper  buckling  load  for  the  additional  load  A P  is 
found  by  setting  d(AP)/db  =  0.  This  yields  the  following  quadratic 
equation  for  the  total  deflection  6^  at  buckling 
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3 &u  "  6a&u  "  bl  +  2a6o  +  2&2  *  ^3/AE)[  (k-(2TQ/L)]  =  0 

The  upper  buckling  load  condition  will  be  reached  by  the  system 
with  no  additional  load  being  required  (AP  =  0)  when  =  bQ,  or 
when 

T0  =  (kL/2)  +  (AE/L2)(6q  -  a)2  (112) 

This  has  an  absolute  minimum  at  TQ  *  kL/2,  The  minimum 
value  that  the  force  in  the  arms  could  have  at  the  upper  buckling 
point  for  any  arbitrary  loading  path  is  therefore  the  same  as  at 
the  energy  load  for  the  fully  elastic  case. 

The  transitional  buckling  condition  is  reached  when  the 
force  in  the  arms  equals  the  Euler  buckling  load,  T  =  T,?.  For  snap 
buckling  to  occur  for  additional  loading  above  P  ,  the  additional 
_oad  must  be  greater  than  AP0,  (110),  The  total  force  in  the  arms 
must  therefore  exceed  kl/2  if  buckling  of  the  arms  is  to  lead  to 
snap  buckling  and  not  to  only  a  bifurcation  of  equilibrium  with  a 
continuously  increasing  deflection.  Fig,  21  shows  the  combinations 
of  the  force  in  the  arms  and  the  displacement  that  will  cause  in¬ 
stability,  The  relative  values  of  the  different  criteria  depend 
upon  the  original  geometry  and  the  elastic  constants,  but  not  on 
the  creep  properties.  The  relation  between  TQ  and  bQ  for  com¬ 
pletely  elastic  deformation  is  given  by  (I4.)  and  is  also  plotted  in 
Fig,  21,  The  system  described  by  Fig,  21  would  buckle  under  elastic 
loading  when  the  elastic  curve  crosses  the  curve  for  the  limit 
point  (the  upper  buckling  load)  condition,  point  ft. 

For  the  creep  problem  the  system  is  first  loaded  elastically 
up  to  some  point  A  along  the  elastic  response  curve  0Q,  The 
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system  then  creeps  under  the  constant  load  with  consequent  changes 
In  the  foroe  In  the  arms  and  in  the  deflection.  The  relation  be¬ 
tween  these  two  quantities  during  creep  (path  AB)  depends  upon  the 
particular  creep  law  and  would  result  in  snap  buckling  when  the 
values  agree  with  one  of  the  buckling  criteria  determined  previous¬ 
ly#  e.g,,  point  B  in  Pig,  21,  Although  the  law  for  metallic  creep 
is  nonlinear,  the  principal  features  of  the  problem  can  be  most 
easily  obtained  from  the  study  of  linear  viscoelastic  behavior. 

The  creep  laws  for  the  spring  and  the  arms  of  the  model 
can  be  taken  as 


.  *  <Vk)  +  Xpo 

(113) 

(1/A)[(T0/E)  +  (lT0) 

(114) 

where  \  and  p,  are  material  constants. 

Upon  differentiating  (3)  (with  the  sign  changed)  with 
respect  to  the  time,  the  strain  rate  In  the  arms  can  be  obtained 
as  a  function  of  the  displacement, 

z  *  (1/L2)b0(a-6C) 
which,  upon  substitution  Into  (111*.),  gives 

tc(a-&0)  *  (l2/a)[(tc/e)  +pTe]  (115) 

Differentiation  of  the  equilibrium  equation  (108)  with 
respect  to  the  displacement  gives 
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dP  dT 

ar  -  <2/L)tV  <**c>  vh 

C  0 


(116) 


since  the  applied  load  PQ  is  oonstant  during  the  creep  process* 


Using  the  relationship 


pc  «  6 


• 

e 

o 


(117) 


in  (113)*  And  solving  for  the  displacement  rate  gives 


dP 


60  -  XPe/[l-(l/k) 


(118) 


dP. 


Substituting  for  PQ  and  «  from  (108)  and  (116)  leads  to  the  fol- 

o 

lowing  expression  for  the  deflection  rate: 


b 


XPe-(2\/L)(a-60)T0 
1-  (2At)[T0-(.-60)^£j 


(119) 


Using  the  relation 


.  .  dT 

T  *  6  nr  - 

o  c  ®T 
o 


equation  (115)  can  also  be  solved  for  the  deflection  rate  to  give 

(120) 


b 


(lV*)t« 


(a-b  )-(L2/AB)v& 

c 

Buckling  will  take  plaoe  when  bQ  becomes  infinite  or  when 
the  denominators  of  (119)  or  (120)  beoomes  zero.  The  buckling 
criterion  is  therefore 

dT, 


l-(2/kLHT0.  (a-60)  ^2)  -  0 

0 

o  dT_ 

(a-6  )  -  (L  / AE)  -jr—  =  0 

o 


(121) 


(122) 
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Solving  (121)  and  (122)  simultaneously  for  Tc  gives  the 
same  instability  locus  in  the  T  f  *>0,  plane  as  was  obtained  pre¬ 
viously  (112)  by  the  additional  loading  approach.  The  equation 
for  the  creep  locus  in  the  T  ,  &c  plane  (path  AB)  can  be  obtained 
by  equating  the  two  creep  rates  in  (119)  and  (120),  This  gives 

dTc  (X/^)(a-&0)[P0-(2Tc/L)(a-60)]^(2TeL/Ak)[T0-(kL/2))  (123) 

^  "  (X/h.)  ( L2/AE )  [  PQ-  ( 2TC/L )  ( a-6 0 )  ]  +(2T0L/Ak )  ( a-&c ) 

The  required  locus  can  be  obtained  by  integrating  (123), 

The  initial  condition  is  the  solution  for  the  initial  elastic  load¬ 
ing  (point  A  of  *'ig.  21),  When  noint  B  is  reached  the  system  will 
snap  through  to  some  point  C  and  creep  again  along  CD,  The  time 
necessary  to  reaoh  B,  the  creep  buckling  time,  can  be  computed 
from  (123)  and  (119)  or  (120), 

The  slope  of  the  creep  locus  can  be  seen  to  be  positive 

over  most  of  the  T  ,  plane.  The  term 

O'  o 

P0-(2T0/L)(a-6c) 

in  (123)  is  equal  to  F0  from  the  equilibrium  equation  (108)  and 

is  therefore  always  positive.  The  quantity  (a-&c)  is  also  positive 

for  6.  <  a.  The  denominator  of  (123)  is  therefore  always  positive 
c 

in  the  region  0  <  <  a  for  all  T„.  The  first  term  in  the 

numerator  of  (123)  is  also  positive  in  this  region.  The  second 
term  is  positive  for  T  >  kL/2  for  all  6  and  negative  for 
Tc  <  kL/2,  The  slope  of  the  creep  locus  in  the  Tc,  &,  plane  will 
therefore  be  zero  or  negative  for  <  kL/2  when 
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(X/|x)(ft-6e)tP0-(2T0/L)(a-6e)]^(?TcL/Ak)[Tc-(kL/2)]  £0  (124) 

or  in  the  region  below  EP  in  Pig.  21,  The  lowest  initial  point 
that  will  lead  to  buckling  will  be  some  point  G  for  which  the 
oreep  locus  will  be  OP  which  passes  through  P  with  a  horizontal 
elope.  Initial  loading  to  any  point  above  G  will  lead  to  buckling. 
Since  Tq  for  point  G  would,  in  general,  be  less  than  the  force  in 
the  arms  at  the  energy  load,  loads  below  the  energy  load  could 
result  in  snap  buckling  of  the  system, 

T,  H.  H.  Plan  and  C.  Y.  Chow  [191  studied  the  creep  buok- 
ling  of  a  simply  supported  low  aroh.  To  simplify  the  analysis 
they  assumed  that  the  stress  distribution  across  the  arch  is  linear, 
and  that  the  stresses  and  displacements  varied  along  the  arch  in 
the  same  manner  as  in  the  elastic  solution.  The  second  of  these 
conditions  is  equivalent  to  assuming  that  the  deformed  shape  is  the 
same  in  the  elastlo  and  creep  preoesses.They  solved  the  problem 
using  one  symmetrical  and  one  asymaetrical  term  in  the  deflection 
function,  Por  symmetrical  oreep  deformations  they  found  the  arch 
to  behave  in  exactly  the  same  manner  as  the  model.  They  used  a 
creep  law  of  the  form 

e  *  ( d/E)  + 

with  n  =  1  and  3. 

As  no  experimental  results  were  obtainable  in  the  litera¬ 
ture, it  was  felt  that  some  experiments  would  be  very  desirable  in 
view  of  the  many  simplifying  assumptions  mc.de  in  the  analyses,  A 
series  of  tests  were  therefore  performed  on  clamped  arches  in  the 
testing  setup  described  previously.  The  material  used  was  5002-0 


aluminum.  This  was  selected  because  at  a  constant  elevated  tem¬ 
perature  its  mechanical  properties  are  fairly  constant  with  time. 
The  specimen  was  mounted  in  the  oven  at  room  temperature.  One 
end  was  left  loose  and  free  to  expand  so  that  no  thermal  stresses 
would  be  induced  upon  heating.  After  the  oven  and  the  specimen 
had  reached  the  test  temperature.  5>00°F,  and  held  there  for 
2-3  hours  to  ensure  uniform  temperature,  the  specimen  was  clamped. 
This  was  made  possible  by  putting  extensions  on  the  clamping  bolts 
reaching  out  of  the  oven.  The  dial  gauges  were  then  lowered  and 
the  dead  load  applied  at  the  arch  center.  In  order  to  read  all 
the  gauges  and  the  time  simultaneously,  the  garages  and  a  clock  were 
all  photographed  at  regular  intervals.  This  method  proved  success¬ 
ful  and  the  deflected  shape  right  up  to  the  time  of  buckling  was 
obtained  for  all  the  tests.  In 

Table  II  the  loads  and  the  geometry  of  the  arches  are  shown  with 
the  time  required  for  buckling  to  occur.  The  geometries  fall  in 
one  group  around  X  *  9  and  one  around  Is  16,  All  the  tests  were 
performed  at  loads  below  the  energy  load  based  on  the  original 
geometry. 

It  is  Interesting  to  note  tests  number  1  and  if  where  the 
geometry  and  the  load  is  the  same  but  the  buckling  times  differ 
by  a  factor  of  ten.  One  of  the  reasons  for  this  can  be  seen  in 
Pig,  22,  where  the  deflection  Just  prior  to  buckling  for  the  two 
tests  are  plotted.  Specimen  number  If,  the  one  with  the  longest 
buckling  time,  is  almost  symmetrical  whereas  number  1  has  a  large 
unsymmetrical  component. 
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Pig,  23  shows  the  deflections  for  specimen  number  7  at  three 
different  times.  In  both  Pigs,  22  and  23  the  importance  of  includ¬ 
ing  a  second  symmetrical  term  in  the  displacement  function  is 
brought  out.  This  is  the  same  as  was  found  for  the  completely 
elastic  case.  The  Inclusion  of  such  terms  would, however, make  the 
creep  analyses  quite  complex,  A  simplification  may  be  obtained  by 
replacing  the  arch  with  a  structure  of  the  kind  shown  in  Pig,  2I4. 
with  six  elastic-creep  hinges.  The  deflected  shape  of  such  a 
system  would  closely  resemble  the  real  shapes  shown  in  Pigs,  22 
and  23  for  symmetrical  deformations.  Pig,  25  shows  a  buckled  arch 
which  shows  clearly  the  sharp  hinges  formed  at  the  center  and  at 
the  ends.  The  possibility  of  an  analyses  on  these  lines  will  be 
investigated  in  a  later  paper. 
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Table  I 

Dimensions  and  Experimental 
Buckling  Loads  of  Clamped  Arches 


Height (inohes) 

r— — "  — . . 

X  «  (p2R/t) 

.P  (lb.) 

0.346 

3.69 

no  buckling 

0.346 

3.69 

9.05 

0.502 

5.35 

12.15 

0.552 

5.89 

13.70 

0.705 

7.52 

18.10 

0.751 

8.01 

20.85 

0.665 

9.44 

21.30 

0.890 

9.50 

23.90 

1.080 

11.51 

28.45 

1.090 

11,62 

29.65 

1.236 

13.15 

36.20 

1,288 

13.75 

33.30 

1.430 

15.25 

42.00 

1.522 

16.25 

39.50 

Horizontal  distance  between  ends  *=  34  Inches 
Width  *  1  inch;  Thickness  ®  3/16  inch 
Material:  2024-T4  Aluminum  alloy 


TABLE  II 


a 


o* 

0 


* 


K 

vr\ 


L 

m  d 


■  H  • 

O  H  It 

B"a» 

SSg 

»  «  Js  •• 
'  -  dfc 
o 


ipl^o 

•°  ®  §8 

9  o  ® 

-  S  A  & 

<rt  OVA  *3 

tj  do  « 

«rt\A  d 
*ca  r-*  ••  8. 

"S  ii  «  « 

O  «r<  E-* 
*tp  «4> 

d  *d  -p  « 

j  • 

Kl**H 


65 


FIG.  I  TYPICAL  LOAD  DEFLECTION  CURVE  FOR  SYSTEM 
THAT  EXHIBITS  SNAP  BUCKLING. 


FIG. 2  MODEL  TO  EXIBIT  SNAP  BUCKLING. 


FIG.  3  LOAD  DEFLECTION  CURVE  FOR  MODEL. 


FIG.  4  GRAPHICAL  DETERMINATION  OF  ENERGY  LOAD  FOR  THE  MODEL. 


FIG. 6  LOAD  DEFLECTION  CURVE  FOR  THE  MODEL 
SHOWING  NON- SYMMETRICAL  BUCKLING  PATH. 


FIG.  II  LOAD  DEFLECTION  CURVE  FOR  ARCH  SHOWING 
ANTI  -  SYMMETRICAL  TRANSITIONAL  MODE. 


SYMMETRICAL  BUCKLING,  WITH  TRANSITIONAL  ANTI -SYMMETRICAL 

TRANSITIONAL  MODE  PRESENT  BUCKLING 

IN  UNSTABLE  REGION 


O2  o 

FIG.  12  VARIATION  OF  BUCKLING  MODES  WITH  X  (X=--^  ) 


FIG.  13  THEORETICAL  AND  EXPERIMENTAL  BUCKLING  VALUES  FOR  CLAMPED  ARCH  SUBJECTED  TO  CENTRAL  LOAD 


FIG.  14  ELASTIC  AND  CREEP  BUCKLING  APPARATUS. 


FIG.  15  APPARATUS  FOR  ELAStlC  AND  CREEP 
BUCKLING  TESTS  OF  ARCHES. 


FIG.  16  INSIDE  OF  OVEN  SHOWING  AUXILIARY  FRAME  AND  LOADED  SPECIMEN. 
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FIG.  17, A  LOAD  DEFLECTION  CURVES  FOR  THEORY  AND  EXPERIMENT,  A=  11.62 


(PR/Et*f/j)XIO 


FIG.  17  B  EXPERIMENTAL  LOAD  DEFLECTION  CURVE 


FIG.  17  C 
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FIG.  18  DEFLECTED  SHAPE  AT  THE  BUCKLING  LOAD. 
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FIG.  19  DEFLECTED  SHAPE  AT  THE  UPPER  AND  LOWER  BUCKLING 
LOAD,  AND  AT  THE  MID- POSITION. 


FIG. 20  LOAD  DEFLECTION  RELATION  FOR  ELASTIC  LOADING 
BEYOND  AN  INITIAL  STRESS  STATE. 


FIG.  24  SIMPLIFIED  ARCH . 


FIG.  25  BUCKLED  CREEP  SPECIMEN  NO.  5  AND  FIXING  BLOCKS  SHOWING 
SHARP  CURVATURE  CHANGES  AT  THE  CENTER  AND  THE  ENDS. 
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Attn:  Technical  Library  U) 

Chief 

Armed  Forces  Special  Weapons 
Project 
The  Pentagon 
Washington  25,  D.C. 

Attn:  Technical  Information 

Division  (2) 

Weapons  Effects  Div.  (1) 
Special  Field  Projects(l) 
Blast  and  Shock  Branch(l) 

Office  of  the  Secretary  of 
the  Army 
The  Pentagon 
Washington  25,  D.C. 

Attn:  Army  Library  11) 

Chief  of  Staff 
Department  of  the  Army 
Washington  25,  D.C. 

Attn:  Development  Branch 

Res.  and  Dev.  Division(l) 
Research  Branch 
Res.  and  Dev.  Division(l) 
Special  Weapons  Branch 
Res.  and  Dev.  Division(l) 

Commanding  Officer 

Engineer  Research  Development 

Laboratory 

Fort  Bel  voir,  Virginia  (1) 
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(2) 


Office  of  the  Chief  of  Ordnance 
Department  of  the  Army 
Washington  25 *  D.C. 

Attn:  Research  and  Materials 
Branch 

(Ord.  Res.  and  Dev.  Div. )  (1) 

Office  of  the  Chief  of  Engineers 
Department  of  the  Army 
Washington  25*  D.C. 

Attn:  ENG-HL  Lib.  Br. ,  Adm.  Ser. 

Div.  (1) 

ENG-WE  Eng.  Div.,  Civil 

Works  (1) 

ENG-EB  Prot.  Constr.  Br., 

Eng*  Div.,  Mil. 

Constr.  (1) 

ENG-WD  Planning  Div.  Civil 
Works  (1) 

ENG-EA  Struc.  Br.,  Eng.  Div., 
Mil.  Constr.  (1) 
ENG-NB  Special  Engr.  Br., 

Eng.  Res.  and  Dev. 
Division  (1) 


Commanding  Officer 
Squier  Signal  Laboratory 
Port  Monmouth,  Hew  Jersey 
Attn:  Components  and  Materials 
Branch  ( 1 ) 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington  25* 

Attn:  Op  37 

Commandant,  Marine  Corps 
Headquarters,  U.S.  Marine  Corps 
Washington  25#  D.C.  (1) 

Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25#  D.C. 

Attn:  Code  312  (2) 

Code  376  (1) 

Code  377  (1) 

Code  U20  (1) 

Code  1+23  (2) 

Code  1*1*2  (2) 


Office  of  the  Chief  Signal  Officer 
Department  of  the  Army 
Washington  25*  D.C. 

Attn:  Engineering  and  Technical 
Division  (1) 

Commanding  Officer 

Watertown  Arsenal 

Watertown,  Massachusetts 

Attn:  Laboratory  Division  (1) 


Chief,  Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25#  D»C. 

Attn:  AE-1* 

AV-34 

AD 

AD-2 

TD-1*2 

RS-7 

rs-B 


(1) 

(1) 

(1) 

(1) 

(1) 

(1) 

(1) 


Commanding  Officer 
Prankford  Arsenal 
Bridesburg  Station 
Philadelphia  37*  Pennsylvania 
Attn:  Laboratory  Division  (1) 

Office  of  Ordnance  Research 

2127  Myrtle  Drive 

Duke  Station 

Durham,  North  Carolina 

Attn:  Division  of  Engineering 

Sciences  (1) 


Chief,  Bureau  of  Ordnance 


Department  of  the  Navy 
Washington  25*  D.C. 

(1) 

Attn:  Ad3 

Re 

(1) 

Res 

(1) 

Reu 

(1) 

ReS5 

(1) 

ReSl 

(1) 

Ren 

(1) 
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thief,  Bureau  of  Yard*  and  Dock* 
•apartment  of  the  Navy 
Washington  25,  D.C. 


Code  D-202 

(1) 

Code  D-202. 3 

(l! 

Code  D-220 

(1) 

Code  D-222 

(1) 

Code  D-410C 

(1) 

Code  D-qi^O 

(1) 

Code  D-500 

(1) 

Commanding  Officer  and  Director 
David  Taylor  Model  Baain 
Washington  7,  D.C. 


Attn:  Code  lij.0 

(1) 

Code  600 

(1) 

Code  700 

(1) 

Code  720 

(1) 

Code  725 

(1) 

Code  731 

(1) 

Code  7^0 

(1) 

Commander 

U.  3.  Naval  Ordnance  Laboratory 
Khite  Oak,  Maryland 
Attn:  Technical  Library  (2) 

Technical  Evaluation 
Department  (l) 

Director 

Materials  Laboratory 

hew  York  Naval  Shipyard 

Brooklyn  1,  New  York  (1) 

Coamanding  Officer  and  Director 
U.  S,  Naval  Electronics  Laboratory 
San  Diego  $2,  California  (1) 

Of flcer-ln -Charge 
Naval  Civil  Engineering  Research 
and  Evaluation  Laboratory 
U.  S.  Naval  Construction 
Battalion  Center 

Port  Hueneme,  California  (2) 

Director 

Naval  Air  Experimental  Station 
Naval  Air  Material  Center 

Naval  Baae 

Philadelphia  12,  Pennsylvania 
Attn:  Materials  Laboratory  (1) 

Structures  Laboratory  Cl) 


(3) 


Officer-in-Charge 
Underwater  Jbq>loslon  Researoh 
Division 

Norfolk  Naval  Shipyard 

Portsmouth,  Virginia 

Attn:  Dr.  A.  H.  Kell  (2) 

Commander 

U,  S.  Naval  Proving  Grounds 
Dahlgren,  Virginia  (1) 

Sup  er  i  nt  enden  t 

Naval  Sun  Factory 

Washington  25,  D.C.  (X) 

Commander 

Naval  Ordnance  Test  Station 
Inyokern,  China  Lake,  California 
Attn:  Physics  Division  (1) 

Mechanics  Branch  (1) 

Commanding  Officer  and  Director 
Naval  Engineering  Experimental 
Station 

Arvupolls,  Maryland  (1) 

Coamanding  Officer 
USNNOEU 

Kirtland  Air  Porce  Base 
Albuquerque,  New  Mexico 
Attn:  Code  20 

XDr.  J.  N.  Bremnen)  (1) 

Superintendent 

Naval  Post  Graduate  School 

Monterey,  California  (1) 

Commandant 
Marine  Corps  School 
Quantico,  Virginia 
Attn:  Director,  Marine  Corps 
Development  Center  (1) 

Coamanding  General 
U.  S.  Air  Porce 
Washington  25,  D.C. 

Attn:  Research  and  Development 
Division  (1) 
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U) 


Commander 

Air  Material  Command 
Wright-Pat ter son  Air  Force  Base 
Dayton,  Ohio 

Attn:  MCREX-B  (1) 

Structures  Division  (1) 

Commander 

U.  S,  Air  Force  Institute  of 
Technology 

Wright-Patterson  Air  Force  Base 
Dayton,  Ohio 

Attn:  Chief,  Applied  Mechanics 

Group  (l) 

Director  of  Intelligence 
Headquarters,  U.  S,  Air  Force 
Washington  25,  D.C, 

Attn:  P.  V.  Branch 

(Air  Targets  Division)  (1) 

Commander 

Air  Force  Office  of  Scientific 
Research 

Washington  25,  D.C. 

Attn:  Mechanics  Division  (1) 

U.  S,  Atomic  Energy  Commission 
Washington  25,  D.C. 

Attn:  Director  of  Research  (2) 

Director 

National  Bureau  of  Standards 
Washington  25,  D.C. 

Attn:  Division  of  Mechanics  (1) 
Engineering  Mechanics 
Section  (1) 

Aircraft  Structures  (1) 

Commandant 

U.S.  Coast  Guard 

1300  E.  Street,  N.  W. 

Washington  25,  D.C, 

Attn:  Chief,  Testing  and 

Development  Division  (1) 

U.  S.  Maritime  Administration 
General  Administration  Office 
Building 

Washington  25,  D.C. 

Attn:  Chief,  Division  of  Pre¬ 
liminary  Design  (1) 


National  Aeronautics  and 
Space  Administration 
1515  H  Street,  N.  w. 

Washington  25,  D.C. 

Attn:  Loads  and  Structures 

Division  (2) 

Director 

Langley  Aeronautical  Laboratory 

Langley  Field,  Virginia 

Attn:  Structures  Division  (2) 

Director 

Forest  Pro  ducts  Laboratory 
Madison,  Wisconsin  (1) 

Civil  Aeronautics  Administration 
Department  of  Commerce 
Washington  25,  D.C. 

Attn:  Chief,  Aircraft  Engineer¬ 
ing  Division  (1) 

Chief,  Airframe  and 
Equipment  Branch  (1) 

National  Sciences  Foundation 
1951  Constitution  Avenue,  N.W. 
Washington,  D.C. 

Attn:  Engineering  Sciences 

Division  (1) 

National  Academy  of  Sciences 
2101  Constitution  Avenue 
Washington  25,  D.C. 

Attn:  Technical  Director, 
Committee  on  Ships 
Structural  Design  (1) 

Executive  Secretary, 
Committee  on  Undersea 
Warfare  (1) 

Professor  Lynn  S.  Beedle 
Fritz  Engineering  Laboratory 
Lehigh  University 
Bethlehem,  Pennsylvania  (1) 

Professor  R.  L.  Bisplinghoff 
Dept,  of  Aeronautical  Engineering 
Massachusetts  Institute  of 
Technology 

Cambridge  39,  Massachusetts  (1) 
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Professor  H.  H,  Bleich 
Department  of  Civil  Engineering 
Columbia  University 
New  York  27,  New  York  (1) 

Professor  B.  A.  Boley 
Department  of  Civil  Engineering 
Columbia  University 
New  York  27,  New  York  (1) 

Professor  G.  P.  Carrier 

Pierce  Hall 

Harvard  University 

Cambridge  38,  Massachusetts  (1) 

Professor  Herbert  Deresiewicz 

Dept,  of  Mechanical  Engineering 

Columbia  University 

632  W.  125th  Street 

New  York  27,  New  York  (1) 

Professor  D.  C.  Drucker 
Division  of  Engineering 
Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  A.  C.  Eringen 
Division  of  Engineering  Sciences 
Purdue  University 
Lafayette,  Indiana  (1) 

Professor  W.  Fltigge 

Dept,  of  Aeronautical  Engineering 

Stanford  University 

Stanford,  California  (1) 

Professor  J.  N,  Goodler 

Dept,  of  Engineering  Mechanics 

Stanford  University 

Stanford,  California  (1) 

Professor  L.  E.  Goodman 
Engineering  Experiment  Station 
University  of  Minnesota 
Minneapolis,  Minnesota  (1) 

Professor  M.  Hetenyi 
The  Technological  Institute 
Northwestern  University 
Evanston,  Illinois 


(5) 


Professor  P.  G.  Hodge 
Department  of  Mechanics 
Illinois  Institute  of  Technology 
Chicago  16,  Illinois  (1) 

Professor  N.  J.  Hoff 

Dept,  of  Aeronautical  Engineering 

Stanford  University 

Stanford,  California  (1) 

Professor  V..  H.  Hoppmann,  II 
Department  of  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  (1) 

Professor  Bruce  G.  Johnston 

University  of  Michigan 

Ann  Arbor,  Michigan  (1) 

Professor  J.  Kempner 
Dept,  of  Aeronautical  Engineering 
and  Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
99  Livingston  Street 
Brooklyn  2,  New  York  (1) 

Professor  H.  L.  Langhaar 
Dept,  of  Theoretical  and 
Applied  Mechanics 
University  of  Illinois 
Urbana,  Illinois  (l) 

Professor  B.  J.  Lazan,  Director 

Engineering  Experiment  Station 
University  of  Minnesota 
Minneapolis  II4.,  Minnesota  (l) 

Professor  E.  H.  Lee 

Division  of  Applied  Mathematics 

Brown  University 

Providence  12,  Rhode  Island  (1) 


(1) 
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Professor  Paul  Lieber 
Geology  Department 
University  of  California 
Berkeley  4,  California  (1) 

Professor  Hsu  Lo 

School  of  Engineering 

Purdue  University 

Lafayette,  Indiana  (1) 

Professor  R.  D.  Mindlin 

Department  of  Civil  Engineering 

Columbia  University 

632  V..  125th  Street 

New  York  27,  New  York  (1) 

Dr.  A.  Nadai 

136  Cherry  Valley  Road 

Pittsburgh  21,  Pennsylvania  (1) 

Professor  Paul  M.  Naghdi 
Mech.  Bngin.,  Mechanics  &  Design 
University  of  California 
Berkeley  4,  California  (1) 

Professor  William  A.  Nash 
Dept. of  Engineering  Mechanics 
University  of  Florida 
Gainesville,  Florida  (1) 

Professor  N.  M.  Newmark,  Head 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  Aria  Phillip* 

Department  of  Civil  Engineering 
15  Prospect  Street 
Yale  University 

New  Haven,  Connecticut  (l) 

Professor  W.  Prager 
Computing  Center 
Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  E.  Reissner 
Department  of  Mathematics 
Massachusetts  Institute  of 
Technology 

Cambridge  39,  Massachusetts  Cl) 


(6) 


Professor  M.  A.  Sadows^ty-^ 
Department  of  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  (1) 

Professor  J.  Stallmeyer 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  Eli  Sternberg 
Division  of  Applied  Mathematics 
Brown  University 
Providence  12,  Rhode  Island  (1) 

Professor  S.  P.  Timoshenko 
School  of  Engineering 
Stanford  University 
Stanford,  California  (1) 

Professor  A.  S.  Velestos 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  Enrico  Volterra 
Dept,  of  Engineering  Mechanics 
University  of  Texas 
Austin,  Texas  (l) 

Professor  Dana  Young 
Yale  University 

New  Haven,  Connecticut  (1) 

Professor  Bernard  Vi.  Shaffer 
Dept,  of  Mechanical  Engineering 
New  York  University 
New  York  53»  New  York  (1) 

Dr.  John  F.  Brahtz 

Southern  California  Laboratories 

Stanford  Research  Institute 

820  Mission  Street  (l) 

South  Pasadena,  California 

Mr.  Martin  Goland,  President 

Southwest  Research  Institute 

8500  Culebra  Road 

San  Antonio,  Texas  (1) 

Mr.  S.  Levy 

General  Electric  Research 
Laboratory 
6901  Elmwood  Avenue 
Philadelphia  42,  Pennsylvania  (1) 


Nonr-562(20)  Distribution  List 


Professor  B.  Budiansky 

Dept,  of  Mechanical  Engineering 

School  of  Applied  Sciences 

Harvard  University 

Cambridge  39*  Massachusetts  (1) 

Professor  George  Herrmann 
Department  of  Civil  Engineering 
Columbia  University 
New  York  27,  New  York  (1) 

Professor  E.  Orowan 
Dept,  of  Mechanical  Engineering 
Massachusetts  Institute  of 
Technology 

Cambridge  39,  Massachusetts  (1) 

Professor  J.  Ericksen 
Mechanical  Engineering  Department 
Johns  Hopkins  University 
Baltimore  18,  Maryland  (1) 

Professor  T.  Y.  Thomas 
Graduate  Institute  for  Mathematics 
and  Mechanics 
Indiana  University 
Bloomington,  Indiana  (1) 

Professor  Joseph  Marin,  Head 
Department  of  Engineering  Mechanics 
College  of  Engineering  and 
architecture 

The  Pennsylvania  State  University 
University  Park,  Pennsylvania  (1) 

Professor  Robert  L.  Ketter 
Department  of  Civil  Engineering 
University  of  Buffalo 
Buffalo  14,  New  York  (l) 

Professor  W.  A.  Shaw 
Department  of  Mechanical  Engineering 
Alabama  Polytechnic  Institute 
Auburn,  Alabama  (1) 

Mr.  K.  H.  Koopman,  Secretary 
Welding  Research  Council  of 
The  Engineering  Foundation 
29  W  39th  Street 

New  York  18,  New  York  (2) 


(7) 


Professor  Walter  T.  Daniels 
School  of  Engineering  and 
Architecture 
Howard  University 
Washington  1,  D.C.  (1) 

Professor  P.  S.  Symonds,  Chairman 
Division  of  Engineering 
Brown  University 
Providence  12,  Rhode  Island  (1) 

Professor  Nicholas  Per rone 
Engineering  Science  Department 
Pratt  Institute 

Brooklyn  5,  New  York  (1) 

Commander 

Wright  Air  Development  Center 
Wright-Patterson  Air  Force  Base 
Dayton,  Ohio 

Attn:  Dynamics  Branch  (1) 

Aircraft  Laboratory  (1) 

WCLSY  (1) 

Dr.  Edward  Wenk,  Jr. 

Executive  Secretary 

Federal  Council  for  Science  &  Technology 
The  White  House 

Washington,  D.  C.  (1) 

Dr.  H.  G.  Hopkins 
War  Office 

Armament  Research  and  Development 
Establishment 
Fort  Halstead 

Sevenoaks,  Kent,  England  (1) 

Professor  J.  E.  Cermak 
Department  of  Civil  Engineering 
Colorado  State  University 
Fort  Collins,  Colorado  (1) 

Professor  Vi.  J.  Hall 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  R.  Muki 

Division  of  Mechanical  Engineering 

Keio  University 

Koganei-shi 

Tokyo,  Japan  (1) 


